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^^ Abstract. In this paper we examine in detail the procedure of averaging of the local field- 



theoretic Poisson brackets proposed by B.A. Dubrovin and S.P. Novikov for the method of 



r\i Whitham. The main attention is paid to the questions of justification and the conditions of 

- -. applicability of the Dubrovin-Novikov procedure. Separate consideration is given to special 

^ features of single-phase and multiphase cases. In particular, one of the main results is the 

^\ insensitivity of the procedure of bracket averaging to the appearance of "resonances" which 

can arise in the multi-phase situation. 
■^ 
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1 Introduction 



As is well-known, the Whitham method [48-50] is associated with slow modulations of periodic 
h>. or quasiperiodic rn--phase solutions of nonlinear systems 

(N 

en F*((^,(^t,(^,,...)=0, i = l,...,n, (p=(^\...,v,"), (1.1) 

which are usually represented in the form 

m 

^ (/p^(2;,t) = $*(k(U)x + c^(U)i + 6>o,U). (1.2) 



Let us note that we will consider here systems with one spatial variable x and one time 

variable t. In these notations the functions k(U) and tj(U) play the role of the "wave numbers" 

and "frequencies" of ?Ti--phase solutions, while the parameters 6q represent the "initial phase 

Cd shifts". The parameters U = {U^ , . . . , U^) can be chosen in an arbitrary way, we just assume 

that they do not change under shifts of the initial phases of solutions Oq. 

The functions <I>*(0) satisfy the system 

F*($,a;"*e.,A:'5*e,5,...) =0, i = l,...,n, (1.3) 

and we have to choose for each value of U some function $(0, U) as having "zero initial phase 
shift". The corresponding set of m-phase solutions of (1.1) can be then represented in the 
form (1.2). For m-phase solutions of (1.1) we have in this case k(U) = (A;^(U), . . . , /c™(U)), 
u;(U) = (a;i(U), . . . , w'"(U)), Oq = {dl, . . .,9}^), where U = (C/\ ...,U^) are the parameters 
of a solution. We will require also that all the functions <I>*(0,U) are 27r-periodic with respect 
to each 0", a = 1, . . . ,m. 



*This paper is a contribution to the Special Issue "Geometrical Methods in Mathematical Physics". The full 
collection is available at http://www.emis.de/journals/SIGMA/GMMP2012.html 
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Consider a set A of functions ^{6 + Oo,\J), depending smoothly on the parameters U and 
satisfying system (1.3) for all U. 

In the Whitham approach the parameters U and Oq become slowly varying functions of x 
and t, U = U(X,T), Oq = 0o(^,T), where X = ex, T = et (e ->■ 0). 

For the construction of the corresponding asymptotic solution the functions U(X, T) must sa- 
tisfy some system of differential equations (the Whitham system). In the simplest case (see [31]), 
we try to find asymptotic solutions 

^\e, X,T) = Y^ ^J,) fS{X^ ^ Q ^^ ^\ ^k (14^ 

fc>0 ^ ^ ^ 

with 27r-periodic in 6 functions ^ (k) satisfying system (1.1), i.e. 
F'((^,e(/3T,evx,..-) = 0, z = l,...,n. 

The function S{X,T) = {S^{X,T), . . .,S'^{X,T)) is called the "modulated phase" of solu- 
tion (1.4). 

Assume now that the function ^(o)(^)^)^) belongs to the family A of m-phase solutions 
of (1.1) for all X and T. We have then 

*(o)(0,x,r) = *(0 + 0o(^,r),u(x,r)), (1.5) 

and 

5?(X,r) = a;°(U(X,r)), Sl(X,r) = A:°(U(X,r)), 

as follows after the substitution of (1.4) into system (1.1). 

In the simplest case the functions ^(^^i^\{9,X,T) are determined from the linear systems 

L]lUA>](^^T)^U((^,X,T) = fi,^{e,X,T), (1.6) 

where L^uj g ]{^:T) is a linear operator defined by the linearization of system (1.3) on the so- 
lution (1.5). The resolvability conditions of systems (1.6) in the space of periodic functions can 
be written as the conditions of orthogonality of the functions f(A;)(^, X, T) to all the "left eigen- 
vectors" (the eigenvectors of the adjoint operator) of the operator L\-yj q AX,T) corresponding 
to zero eigenvalue. 

We should say, however, that the resolvability conditions of systems (1.6) can actually be quite 
complicated in general multi-phase case, since the eigenspaces of the operators -^[u.eo] ^^^ ^fu e 1 
on the space of 27r-periodic functions can be rather nontrivial in the multi-phase situation. Thus, 
even the dimension of the kernels of -C/ru.eo] ^"^^ ^\r\ e 1 ^^^ depend in a highly nontrivial way 
on the values of U. In general, the picture arising in the U-space can be rather complicated. As 
a result, the determination of the next corrections from systems (1.6) is impossible in general 
multiphase situation and the corrections to the main approximation (1.5) have more complicated 
and rather nontrivial form [4-6]. 

These difficulties do not arise commonly in the single-phase situation [m = 1) where the 
behavior of eigenvectors of -^[u.eo] ^^^^ [u 6» ii ^s a rule, is quite regular. The resolvability 
conditions of system (1.6) for k = 1 

^}[u,«o](^'^)*{i)(^'^'^) = /(i)(^'^'^) (1-7) 

with relations kx = oox define in this case the Whitham system for the single-phase solutions 
of (1.1) which plays the central role in considering the slow modulations. 
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For the multi-phase solutions the Whitham system is usually given by the orthogonality 
conditions of the right-hand part of (1.7) to the maximal set of "regular" left eigenvectors 
corresponding to zero eigenvalues which are defined for all values of U and depend smoothly 
on U. As a rule, this system is equivalent to the conditions obtained by the averaging of some 
complete set of conservation laws of system (1.1) having the local form 

Pi'iip, ^t, ^x,-..) = QU^, ^t, <^x, . . . ), U=l,...,N. 

The Whitham system is written usually as a system of hydrodynamic type 

u^ = y;(u)t/^ (1.8) 

and gives the main approximation for the connection between spatial and tiine derivatives of the 
parameters U^{X,T). The variables T and X represent the "slow" variables T = et, X = ex, 
connected with the variables t and x by a small parameter e. Thus, the Whitham system (1.8) 
represents a homogeneous quasi-linear system of hydrodynamic type connecting the derivatives 
of slow modulated parameters. 

As mentioned above, the construction of the asymptotic series (1.4) in the multi-phase case 
is impossible in general situation (see [4-6]). Nevertheless, the Whitham system (1.8) and the 
leading term of the expansion (1.4) play the major role in consideration of modulated solutions 
also in this case, representing the main approximation for corresponding modulated solutions. 
The corrections to the main term have in general more nontrivial form than (1.4), but they also 
tend to zero in the limit e — )• [4-6]. 

Let us give here just some incomplete list of the classical papers devoted to the foundations of 
the Whitham method [1,4-8,14-16,22-25,28,29,31,39,40,42,48-50]. We wiU be interested here 
only in Hamiltonian aspects of the Whitham method. In the remaining part of the Introduction 
we will give the definition of the "regular" Whitham system for the complete regular family of 
m-phase solutions which will be used everywhere below. 

Let us use for simplicity the notation A both for the family of the functions ^{0 + 9q, U) and 
the corresponding family of ?T2-phase solutions of system (1.1), such that we will denote by A both 
the families of the functions ^{0 + 6q, U) in the space of 27r-periodic in all 0" functions (f{0) and 
¥'[U,0o](^) = *(k(U)x + Oq, U). We will assume everywhere below that the family A represents 
a smooth family of m-phase solutions of system (1.1) in the sense discussed above. 

It is generally assumed that the parameters k°', u" are independent on the family A, such that 
the full family of the m-phase solutions of (1.1) depends on A^ = 2m + s (s > 0) parameters U^ 
and m initials phase shifts 9q. In this case it is convenient to represent the parameters U in the 
form U = (k, u, n), where k represents the wave numbers, u: are the frequencies of the m-phase 
solutions and n = (n^, . . . ,n^) are some additional parameters (if any). 

It is easy to see that the functions ^ga(d + 0O)k, cj, n), a = 1, . . . ,m, ^^n^iO + Oq, k, o;, n), 
/ = 1, . . . , s, belong to the kernel of the operator L\^ t^ n e I- -'■^ *^^ regular case it is natural 
to assume that the set of the functions (^e^j ^nO I'cpresents the maximal linearly independent 
set of the kernel vectors of the operator L regularly depending on the parameters (k, a;, n). For 
the construction of the "regular" Whitham system we have to require the following property of 
regularity and completeness of the family of m-phase solutions of system (1.1). 

Definition 1.1. We call a family A a complete regular family of m-phase solutions of sys- 
tem (1.1) if: 

1) the values k = (fc^, . . . , A;'"), u = (o;^, . . . , w"^) represent independent parameters on the 
family A, such that the total set of parameters of the ?Ti--phase solutions can be represented 
in the form (U, Oq) = (k, u, n, Oq); 
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2) the functions ^Qa{0 + 0O)k, a;,n), ^^^{9 + 0O)k, a;,n) are linearly independent and give 
the maximal linearly independent set among the kernel vectors of the operator L\]^ a; n e 1 ' 
smoothly depending on the parameters (k, o^, n) on the whole set of parameters; 

3) the operator i^-rj^.^^ „ q i has exactly m + s linearly independent left eigenvectors with zero 
eigenvalue 

fcfgj {e + 0o) = ^f^?^,„] (0 + 0o), g = 1, . . . , m + s, 

among the vectors smoothly depending on the parameters (k, a^,n) on the whole set of 
parameters. 

By definition, we will call the regular Whitham system for a complete regular family of 
m-phase solutions of (1.1) the conditions of orthogonality of the discrepancy f(i)(0, X, T) to the 

functions '^[u(x,T)] (^ + ^o(^, T)) 

■■■I '^[U(X,T)].(^ + ^o(X,r))/(\)(0,X,r)^^-^=O, q=l,...,m + s, (1.9) 

with the compatibility conditions 

k^ = uj\. (1.10) 

System (1.9), (1.10) gives ra + (m + s) = 2m + s = A^ conditions at every X and T for the 
parameters of the zero approximation *(o)(^)^i^)- 

It is well known that the Whitham system does not include the parameters Oq (X, T) and 
provides restrictions only to the parameters U'^{X,T) of the zero approximation. Let us prove 
here a simple lemma which confirms this property under the conditions formulated above^. 

Lemma 1.1. Under the regularity conditions formulated above the orthogonality conditions (1.9) 
do not contain the functions 6q{X, T) and give constraints only to the functions U^{X, T), having 
the form 

with some functions Ci. (U), Dy {\J). 

Proof. Let us write down the part f^^ of the function fn), which contains the derivatives 
elj,{X,T) and elx{X,T). We have 



/(i)(0, X, T) = -— (*(o), . . . ) ^\,^At - ^ (*(0), • • • ) "^UeAx 
dF 



dF' 






^This simple fact was present in the Whitham approach from the very beginning (see [31,48-50]). In fact, 
under various assumptions it can be also shown that the additional phase shifts 9q {X, T) can be always absorbed 
by the functions S°'{X,T) after a suitable correction of initial data (see, e.g., [9,23,24,35]). It should be noted, 
however, that the corresponding phase shift can play rather important role in the weakly nonlocal case [39] (see 
also [9,33]). 
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,k"',co\...,u;"',n\. 



,n 



Let us choose the parameters U in the form 

V = {k\...,k"',co^ 

We can write then 

d 



fii){0,X,T) 



dujP 



ri^{e + eo,v), 



+ L 



+ 



d 



F'(*(6' + 6>o,U) 



+ L 



d_ 
d 



^dkP 



$J(0 + 6>o,U) 



'OT 



'ox- 



The total derivatives dF'^ /duj'^ and dF^ jdk^ are identicahy equal to zero on A according 
to (1.3). We have then 



^2-K j'2-K (Y'^R 



since I'^\jj(xt)]^^ ~^ ^o{X,T)) are left eigenvectors of L with the zero eigenvalues. 

It is not difficult to see also that all the Oq{X,T) in the arguments of $ and k^'^' disappear 
after the integration with respect to 0, so we get the statement of the lemma. ■ 

We can claim then that the regular Whitham system has the following general form 






dU' 



dU" 



X' 



a = 1, . . . ,m, 



Cl^^Um = D\:^\\J)Uh q=l,...,m + s. 



:i.iii 



Let us note that according to our assumptions we have here rank | |9A;"/9C/^| | = m. In generic 
case the derivatives C/^ can be expressed in terms of U'^ and the Whitham system can be written 
in the form (1.8). 



2 Lagrangian and Hamiltonian formulations 
of the Whitham method 

Together with the formulation of Whitham's method the Lagrangian structure of the equations 
of slow modulations was proposed [48-50]. The method of averaging of Lagrangian function 
introduced by Whitham can be formulated in the following way. We assume that the original 
system (1.1) is lagrangian with the local action of the form 



S 



L{(^,Lpf,'fx-,ftt-,'fxt-,'fz 



] dxdt, 



such that the functions F* have the form 

i 6S _ dL _ d dL _ d dL 

6ip^{x,t) dip^ dtdifl dx difl. 



+ ■■■ ■ 



Let us assume here for simplicity that the parameters (k, uj) = {k^, . . . , k^,uj^, . . . , uj^) give 
the complete set of independent parameters on the family of 77i-phase solutions (excluding the 
initial phase shifts), such that the number of parameters C/^ is equal to 2m. 

The linearized operator L^]^,^ g ]iX,T) in (1.6) is given now by the distribution 



-^}[k,a;,0o] 



{o,e') 



6^S 



S<^i{e)6'^j{6') 
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where 



is a self-adjoint operator. 

Throughout the paper we wih always understand the integration with respect to 6 as the 
averaging procedure. For this reason, all the integrals over d™0 will be defined with the factor 
l/(27r)'". In particular, we will also assume that the variational derivatives of the type 6S/Sip'^(9) 
are defined as 



Jo VW ^' '(2vr)- 

on the space of 27r-periodic in functions. 

We also define here the delta function 5{6 — 6') and its higher derivatives 5gaj^ Qas{6 — 0') 
on the space of 27r-periodic functions by the formula 



/•27r (V^f)' 

■ I (5e^i...ga, {e - 9')^{9') = i;gc^...g^s (6). 

Jo l^^j 



/•27r f2TT A^T^a' 

Jo 



The functions ^ga, a = 1, . . . ,m, represent both the left and the right eigenfunctions of the 
operator L\^^^ ^ JX, T), corresponding to the zero eigenvalue. 

Under the assumption that the family of the m-phase solutions A is a complete regular 
family of 77i-phase solutions of (1.1) we assume that the functions ^ga [9 + 6q, k, u) are linearly 
independent and give the maximal linearly independent set among the kernel vectors of the 
operator -^Vj^.^^ i smoothly depending on the parameters (k, cj). The regular Whitham system 
is given then by the conditions k^ = w^ and m conditions of orthogonality of the function 
f(i)(0,X,r) to the functions ^ec{9 + 9o,k,u). 

According to the Whitham procedure the Whitham system on the parameters (k, a;) is ob- 
tained from the condition of extremality of the action 



^^'HS] = JJJ^ ■■■/o ^(*'^?*^"'4*.''.---)(^^dXdT (2.1) 

under the conditions A;° = 5^, oj°^ = S^. 

The conditions k^ = oj^ and 6T,/6S°'{X,T) = give a system of 2m equations on the 
parameters (k, cj). 

It is not difficult to see that the system given by the variation of the "averaged" action 
coincides with the conditions of orthogonality of the function {n\{9,X,T) to the functions 
^ga{9 + 0O)k, cj). Indeed, let us consider the action 

= SW [S, <p] + eS(i) [S, (^] + e2s(2) [S, v:>] + • • • 
defined on the functions (f{9,X,T), 27r-periodic in each 6". Taking into account the relation 

6SHX,T)-' i ■■■io ^^"^'''^''^V(^,^,T)(2vr)- 
and the invariance of the action with respect to the shifts 
S(X,T)^S(X,r) + AS, 
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it is easy to see that 

f27r |.27r g^{0) ^r. 



JO JO 



6^'{e,X,T){27ry 



'"'"' '"" ''' '^(O.X.T)^^^^^ (2,2) 



«S"(X,T) 7„ io ^" ' ' ■ '«»j'(e,x,r)(2,r) 

etc. 

Substituting the functions (p{0,X,T) in the form (f{e,X,T) = *(0 + 0o, Sx, St) in the 
relations above, we can see that we have to include now the additional dependence of the 
functions if{0,X,T) on Sx and Sy in relation (2.2). However, due to the relation 





6ipHe,x,T) 

on the family A, the relation (2.2) will not change in this situation. Taking also into account 
the equality 

/(\)(0,x,r)-^-^^^^^ 

we get the required statement. 

Under the assumption of the completeness and regularity of the family A we can see then 
that the averaged action (2.1) defines a lagrangian structure of the regular Whitham system 
in general multiphase case. We should say also that the cases with additional parameters n, 
as a rule, can be also included into the scheme described above with the aid of the Whitham 
"pseudo-phases" [49]. Let us note also that different questions connected with the justification 
of the averaging of Lagrangian functions in different orders can be found in [7] . 

Another approach to the construction of the regular Whitham system is connected with the 
method of averaging of conservation laws. According to further consideration of the Hamiltonian 
structure of the Whitham equations we will assume now that system (1.1) is written in an 
evolutionary form 

ipi = F'{(f,ip^,(p^^,...). (2.3) 

The families of the ?TT,-phase solutions of (2.3) are defined then by solutions of the system 

u;''^U=r{^,k^<fe^,...) (2.4) 

on the space of 27r-periodic in each 0" functions (p{0). 

We will assume that the conservation laws of system (2.3) have the form 

Pt (¥', <^x, <^xx, ..■) = Ql{>^, >^x, fxx, ■■■), 
such that the values 

" + 00 



P'{ip,^x,^xx,---)^X 

-oo 



represent translationally invariant conservative quantities for the system (2.3) in the case of 
the rapidly decreasing at infinity functions <^{x). We can also define the conservation laws for 
system (2.3) in the periodic case with a fixed period K 



1 f^ 



A.Ya. Maltsev 



or in the quasiperiodic case 



1 /"^ 

^"^ = J™ ^^ / P"" ('^' <^^' '^^^' • • • ) dx. 



It is natural also to define the variational derivatives of the functionals I'^ with respect to the 
variations of '^{x) having the same periodic or quasiperiodic properties as the original functions. 
Easy to see then that the standard Euler-Lagrange expressions for the variational derivatives 
can be used in this case. 

Let us write the functionals I^ in the general form 

P" {'P,'Px,^xx,---)<^x (2.5) 

assuming the appropriate definition in the corresponding situations. 

Let us define a quasiperiodic function '^{x) with fixed quasiperiods (/c^, . . . , fc™) as a smooth 
periodic function (p(0) on the torus T™ which is restricted on the corresponding straight- line 
winding 



Let us define the functionals 

on the space of 27r-periodic in 6 functions. 
It's not difficult to see that the functions 



(2.6) 



.H ,. . . . sr 



C.[u](^ + ^o)-^^^^ 



(2.7) 



(^(0)=*(6»+0o,U) 



represent left eigenvectors of the operator L^-ijj g ^ with zero eigenvalues regularly depending on 
parameters U on a fixed smooth family A. 

Indeed, the operator L\jj g , is defined in this case by the distribution 



L-,,,,,(e.e')^siu,-6,.(e-e')-'-'"'-^-"^"' 



V(^') 



¥>(0)=#(0+6»o,U) 



We have 

f27r /•27r 






for any translationally invariant integral of (2.3). Taking the variational derivative of this 
relation with respect to ^■'{0') on A we get the required statement. 
Thus, we can write 

q 

with some smooth functions Cg(U) on a complete regular family A. 

For the construction of the regular Whitham system on a complete regular family of m-phase 
solutions of (2.3) we need a sufficient number of the first integrals (2.5) such that the values of 
the functionals J^ on A represent the full set of parameters V^ = J'^\a- Besides that, we should 
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require that the maximal hnearly independent subset of the functions (2.7) give a complete set 
of linearly independent left eigenvectors of the operator L\jjq i with zero eigenvalues among 
the vectors regularly depending on the parameters U on the family A. 

Coming back to the definition of a complete regular family of ?7i-phase solutions of sys- 
tem (2.3) we can see that in the case of a complete regular family A the number of linearly 
independent vectors (2.7) on A is always finite. More precisely, if A^ = 2m + s is the number of 
parameters of m-phase solutions of (2.3) (excluding the initial phase shifts) then for a complete 
regular family of m-phase solutions we require the presence of exactly m + s = N — m left eigen- 
vectors Kr5j, (0 + Oq) with zero eigenvalues, regularly depending on parameters, in accordance 
with the number of the vectors ^e^, $^;. Thus, according to Definition 1.1, we assume here 
that the number of linearly independent vectors defined by formula (2.7) is exactly equal to 
m + s = N — m for a complete regular family A. 

We should note that the conditions on the variational derivatives of J^ formulated above do 
not contradict to the condition that the values J'^ {u = 1, . . . , N) can be chosen as parameters U^ 
on the family of m-phase solutions. Indeed, the definition of J" (2.6) explicitly includes the 
additional m functions A;", which provide the necessary functional independence of the values 
of J^ on A. In other words, we can use the Euler-Lagrange expressions for the variational 
derivatives of I^ only on subspaces with fixed quasiperiods (fc^, . . . , k'^). The variation of the 
quasiperiods gives linearly growing variations which do not allow to use the Euler-Lagrange 
expressions. 

Moreover, under the assumptions formulated above, we can show that the condition of the 
completeness of the variational derivatives (2.7) of the functionals J'^ in the space of regular left 
eigenvectors of the operator L\-yj ^ i with zero eigenvalues follows in fact from the condition that 
the values U" = J^\k can be chosen as the full set of parameters (excluding the initial phase 
shifts) on the family A. 

Let us make the agreement that we will always assume here that the Jacobian of the coordi- 
nate transformation 

(k,Cc.,n)^(C/i,...,f/^) 

is different from zero on A whenever we say that the values [/'^(k,a;,n) represent a complete 
set of parameters on A (excluding the initial phase shifts) . 

Under the conditions formulated above let us prove here the following proposition. 

Proposition 2.1. Let h. he a complete regular family of m-phase solutions of system (2.3). Let 
the values {U^, . . . , U^) of the functionals {J^, . . . , J^) (2.6) give a complete set of parameters 
on A excluding the initial phase shifts. Then: 

1) the set of the vectors 

{*^a(0 + 0o,k,a^,n), *„!(0 + 0o,k,cj,n), a = l,...,m, I = l,...,s} 

is linearly independent on A; 

2) the variation derivatives Q^JjJO + Oq), given by (2.7), generate the full space of the regular 
left eigenvectors of the operator L'^-iu g ^ with zero eigenvalues on the family A. 

Proof. Indeed, we require that the rows given by the derivatives 



/dU^ dU^\ fdU^ dU^ \ 
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are linearly independent on A. Using the expressions 



a = i, ... ,771, 



9n' 






on A, we get that the set {^a;'^) ^n'} i^ linearly independent on A and the number of linearly 
independent variation derivatives (2.7) is not less than m + s. 

We obtain then that the variation derivatives (2.7) generate in this case a space of regular 
left eigenvectors of the operator L\-^J g i with zero eigenvalues of dimension (m + s). ■ 

As a remark, let us note here some general fact associated with multiphase solutions of partial 
differential equations. As is well known, the presence of families of multiphase quasiperiodic 
solutions, as a rule, is connected with integrability of system (1.1) by the inverse scattering 
methods. 

The first TTi-phase solutions for the KdV equation given by the Novikov potentials were 
introduced exactly as the functional families represented by the extremals of linear combinations 
of some set of higher integrals of the system, i.e. the families where the variational derivatives 
of integrals of the set become linearly dependent. 

If we have a natural hierarchy of the first integrals and commuting flows of an integrable sys- 
tem, some of the first integrals of the hierarchy (I^, . . . , I'^) are usually used for the construction 
of m-phase solutions. So, the functions given by the conditions 

ci6I^ + ■■■ + Cg6P = 

for all possible {ci, . . . ,Cq) form a complete family of m-phase solutions of the integrable sys- 
tem [41]. 

Thus, for Novikov potentials we have q = m+2 while the dimensions of the families of m-phase 
solutions for KdV (excluding initial phase shifts) are equal to 2m+l. The first 2m +1 of integrals 
of KdV (/^, . . . , /2m+i^ ^^^ j^g used for the construction of parameters {U^, . . . , f/^^+i) on the 
families of ?n,-phase solutions of KdV. The number of linearly independent variational derivatives 
of these functionals on the family of ?7i-phase solutions is exactly equal to m + 1. It's not difficult 
to show also that the variational derivatives of all the higher integrals of the KdV hierarchy are 
given by linear combinations of the variational derivatives of the set (/^, . . . , I'^) on the families 
of m-phase solutions. 

The construction proposed in [41] in fact is used without substantial changes for many sys- 
tems that are integrable by the inverse scattering methods, and represents the basic scheme 
for constructing of m-phase solutions of integrable systems. This circumstance gives therefore 
a convenient method of checking the above relations for most specific examples. 

Let us prove here the following lemma, which we will need in further considerations. 

Lemma 2.1. Let the values U^ of the functionals J^ on a complete regular family of m-phase 
solutions A he functionally independent and give a complete set of parameters [excluding initial 
phase shifts) on A, such that we have A;" = k'^{U^, . . . , f/^). Then the functionals k°'{J^, . . . , J^) 
have zero variational derivatives on A. 

Proof. As we have seen, the conditions of the lemma imply the existence of m independent 
relations 






= 0, a = l,...,m, (2.9) 



on A. 
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For the corresponding coordinates V^ on A this imphes the relations 

N m 

J^A^(U)dt/^ = 5^^(")(U)dfc'^(U) 

for some matrix /il" (U). 

Since W^ provide coordinates on A the matrix iJ,g(\J) has the fuU rank, and therefore in- 
vertible. We can then write 

m N 



dfc" = E(A"')(«)(U)E^^"^(U)dC/^ 



a=l u=l 

The assertion of the lemma follows then from (2.9). ■ 

The regular Whitham system in the described approach can be written as 

{pnT = {Qnx, u = i,...,N, (2.10) 

where (...) denotes the averaging operation on A defined by the formula 



(/(<^,(p,,...))^y"^••y""/(*,A;/^*,,,... 



{IttY 



Let us prove here the following lemma about the connection between the systems (2.10) 
and (1.11). 

Lemma 2.2. Let the values V^ of the junctionals J^ on a complete regular family of m-phase 
solutions A be functionally independent and give a complete set of parameters on A excluding 
the initial phase shifts. Then the system (2.10) is equivalent to (1.11). 

Proof. Let us introduce the functions 



nr(o(¥''¥'x,..-: 



d^l 



for / > 0. 

Using the expression for the evolution of the densities P'^{(p,eipx, ■ ■ ■) we can write the 
following identities 

Pi'{y^,e^x,...)=Y,e'U^(^i){^,e^x,...){F\y^,e^x,...))i^^eQ'<,i^,e^x,...). (2.11) 

l>0 

To calculate the values e{Q'^)x let us put now 

^i{e,x) = ¥(^^ + e,uix)Y (2.12) 

where 5f = fc"(U(X)). 

The operator ed/dX acting on the functions (2.12) can be naturally represented as a sum of 
k"'d/d9°' and the terms proportional to e. So, any expression f{^, ^fx, • • • ) on the submani- 
fold (2.12) can be naturally represented in the form 

fi<p,e^x,...) = ^e'f[i][^M, 

l>0 
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where /[;] [$, U] are smooth functions of ($, ^ga, ^jj^, . . . ) and (U, Ux, Uxx, ■ ■ ■), polynomial 
in the derivatives (Ux, ^xx, • • • ); ^^"^ having degree I in terms of the total number of derivations 
of U w.r.t. X. Note also that the functions $ appear in /m with the phase shift S(X)/e according 
to (2.12). The common phase shift is not important for the integration with respect to 6, so 
let us assume below that the phase shift S(X)/e is omitted after taking all the differentiations 
with respect to X. 

According to (2.11) and (2.4) we can write 



2-7T f2TT 



eiQlx = e I • • • / Q^[i] ^2^^ 



^ d'^e 



JO 

2lT 1-2-K 

i 
i(/)[l]-''«X[l 

«>0 



JO JO i^n \ 



(27r)™ 



(2^- 

It is not difficult to see also that for arbitrary dependence of parameters U of T, the derivative 
of the average {P'^) w.r.t. T can be written as 

/•2tt r'2iT A^^ft 

^u JO ,>Q V ; 

Now, we can write the relations {P'^)t = {Q")x as 

f>27r r2iT , ^ H'^fl 






l>0 

+ nj(;)[o]^^'''^ • • • k'^'~^^x^ef9^i-e^i-^ + n^(0[o]^^'^^ ' ' ' ^'^'^'^^ %!3e"'i...e'"-ix[i] 

The last three terms in the right-hand part represent the integral of the value 

Y.^^ {^mo]'^l^ix[i]+^mi]'^l^ix[o])=^^dPli^/90^ 

l>0 

and are equal to zero. The remaining terms after integration by parts can be written in the form 

f2n r2TT 



jJ-'-jJi 4u(x)] (^) [mo, U{X)) - F^] {0, X) 



(2^ = °' 



.(u) 



where the values Qxj(x)l(^) ^^^ given by (2.7). 
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Consider the convolution (in ly) of the above expression with the values dk'^/dU'^. The 
expressions 

^(U(A-))cg„„«») 

are identically equal to zero according to Lemma 2.1. 
From the other hand we have 

-.1 ■■■I E^^''---^''-'%.e^....e^^-^^h(^^^^^^^^ 






dU" dk^ 



r[^M 



= ^1, (2.13) 

¥j(0)=*(0,U) 



since the variations of the functions $ are insignificant for the values of k'^ according to 
Lemma 2.1. 

We get then that conditions (2.10) imply the relations k^ = oj^-, which are the first part of 
system (1.11). 

Now the conditions 

r-27r r27r 






express the conditions of orthogonality of the vectors (2.7) to the function —^t + -^fib which 
coincides exactly with the right-hand part of the equation (1.7) in our case. Since the linear 
span of the vectors (2.7) coincides with the linear span of the complete set of the regular left 
eigenvectors of the operator L\jjg AX,T) with zero eigenvalues, we get that system (2.10) is 
equivalent to system (1.11). ■ 

Let us note here that it follows from Lemma 2.2 that systems (2.10), obtained from different 
sets of conservation laws are equivalent to each other. In other words, if system (2.3) has 
additional conservation laws then their averaging gives relations following from system (2.10). 

Let us note also that the justification questions discussed above, as a rule, can be considered 
in a simpler way under additional assumptions about the next corrections to the main approxi- 
mation (1.5) (see, e.g. [7]). We note again that here we don't make any additional assumptions 
of this kind and consider the regular Whitham system as an independent object that is associa- 
ted only with description of the main approximation (1.5). As we have said, we will follow this 
approach everywhere in the paper. 

The Hamiltonian properties of systems (2.10) and more general systems (1.8) play very im- 
portant role in their consideration. The general theory of systems (1.8), which are Hamilto- 
nian with respect to local Poisson brackets of hydrodynamic type (Dubrovin-Novikov brackets) 
was constructed by B.A. Dubrovin and S.P. Novikov. Let us give here a brief description of 
the Dubrovin-Novikov Hamiltonian structures and of the properties of the corresponding sys- 
tems (1.8). 

The Dubrovin-Novikov bracket on the space of fields {U^{X), . . . , U {X)) has the form 

{U''{X),U^'{Y)]=g^^'{V)6'{X-Y) + ^^'{\J)Ul5{X-Y), u,fi = 1, . . . ,N. (2.14) 
The Hamiltonian operator corresponding to (2.14) can be written in the form 
d 



j.M = ^.M(u)-^ + 6-/^(u)c/; 



■7 



As was shown by B.A. Dubrovin and S.P. Novikov [14-16], expression (2.14) with non- 
degenerate tensor g'^>^{\J) defines a Poisson bracket on the space of fields U(X) if and only if: 
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1) tensor g'^^(U) gives a symmetric flat pseudo-Riemannian metric with upper indexes on 
the space of parameters {U^, . . . , U^); 

2) the values 

where g'^ {\J)gxfj_{\J) = 6^^, represent the Christoffel symbols for the corresponding met- 
ric gu^,{'V)■ 

As follows from the statements above, every Dubrovin-Novikov bracket with non-degenerate 
tensor g'^^(U) can be written in the canonical form [14-16] 

{n''{X),n^{Y)} = e''6''^S'{X - Y), e"" = ±1, 

after the transition to the flat coordinates n'^ = n'^(U) for the metric giy^(JJ). 
The functionals 

-00 



/+00 
n''{X)dX 
-00 



represent the annihilators of the Dubrovin-Novikov bracket while the functional 

TV 



/+00 1 J' 
-5]e^K)2(X)dX 



represents the momentum functional for the bracket (2.14). 

The Hamiltonian functions in the theory of brackets (2.14) are represented by the functionals 
of hydrodynamic type, i.e. 



/+00 
/i(U)dX 
-00 



The bracket (2.14) has also two other important forms on the space of U(X). One of them 
is the "Liouville" form [14-16] having the form 

{u-[x), [/^(y)} = (7^'^(u) + 7^'^(u))5'(x - y) + ^uk5{x - v) 

for some functions ^^^^(\J). 

The "Liouville" form of the Dubrovin-Novikov bracket is called also the physical form and 
corresponds to the case when the integrals of coordinates U'^ 



/+00 
U''{X)dX 
-00 



commute with each other. 

Another important form of the Dubrovin-Novikov bracket is the diagonal form. It corre- 
sponds to the case when the coordinates U^ represent the diagonal coordinates for the met- 
ric 5i/^(U) and the tensor g^^{\J) in (2.14) has a diagonal form. This form of the Dubrovin- 
Novikov bracket is closely connected with the integration theory of systems of hydrodynamic 
type 

c/^ = y;(u)[/j^ 
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which can be written in the diagonal form 

[/^ = V''{V)U^ (2.15) 

(no summation) and are Hamiltonian with respect to bracket (2.14). 

It was conjectured by S.P. Novikov that aU the systems of hydrodynamic type having the 
form (2.15) and Hamiltonian with respect to any bracket (2.14) are integrable. This conjecture 
was proved by S.P. Tsarev in [46] where the method (the "generalized hodograph method") of 
integration of these systems was suggested. However, the method of Tsarev proved to be appli- 
cable to a wider class of diagonalizable systems of hydrodynamic type which was called by Tsarev 
"semi-Hamiltonian" . As it turned out later in the class of "semi-Hamiltonian systems" fall also 
the systems Hamiltonian with respect to generalizations of the Dubrovin-Novikov bracket — the 
weakly nonlocal Mokhov-Ferapontov bracket [38] and the Ferapontov bracket [18,21]. Various 
aspects of the weakly nonlocal brackets of hydrodynamic type are discussed in [18-21,36,38,43]. 
Let us note also that the generalization of the Dubrovin-Novikov procedure for weakly nonlocal 
case was proposed in [34] . 

Let us describe now the procedure for constructing the Dubrovin-Novikov bracket for the 
Whitham system in the case when the original system (2.3) is Hamiltonian with respect to 
a local field-theoretic bracket 

Wix),ip^{y)} = Y,B'(^k)(^,V>.,...)5^'Hx-y) (2.16) 

fe>0 

with the local Hamiltonian of the form 

H = J PH{^,^^,...)dx, (2.17) 

which was suggested by B.A. Dubrovin and S.P. Novikov [14-16]. 

Method of B.A. Dubrovin and S.P. Novikov is based on the existence of N (equal to the 
number of parameters V^ of the family of m-phase solutions of (2.3)) local integrals 

r = Jp''iip,^,,...)dx, (2.18) 

which commute with the Hamiltonian (2.17) and with each other 

{I\H} = 0, {F,lf'} = (2.19) 

and can be described as follows. 

We calculate the pairwise Poisson brackets of the densities P'^ having the form 

{P^x),pf^{y)} = ^ Al^i^, (^., . . . )5^''\x - y), (2.20) 

fc>0 

where 

according to (2.19). 

The corresponding Dubrovin-Novikov bracket on the space of functions U(X) has the form 

{U-'iX), U'^iY)} = (^n(U)5'(X -Y) + -^^Ul6{X - Y). (2.21) 
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Let us remind that we assume that the parameters U^ coincide with the values of the func- 
tionals I'^ defined on the corresponding quasiperiodic solutions of the family A 

C/- = (P'^(x)). 

The Whitham system (2.10) is Hamiltonian with respect to the Dubrovin-Novikov brac- 
ket (2.21) with the Hamiltonian 

/+00 
{Ph) (U(A:)) dX (2.22) 

-oo 

The proof of the Jacobi identity for the bracket (2.21) was suggested in [32] under certain 
assumptions about the family of TTi-phase solutions of (2.3). Besides that, it was shown in [37] 
that the Dubrovin-Novikov procedure is compatible with the procedure of averaging of local 
Lagrangian functions when carrying out of both the procedures is possible. 

Let us note also that in [44] all the local brackets (2.14) for the Whitham equations for KdV, 
NLS, and SG equations were found. Besides that, in [2] the hierarchies of the weakly nonlocal 
Hamiltonian structures for the Whitham systems for KdV were represented. 

In this paper we give the most detailed discussion and justification of the Dubrovin-Novikov 
procedure separately for the single-phase and the multiphase cases. In the multiphase case we 
will show also that the justification of the Dubrovin-Novikov procedure can be done under 
weaker assumptions than in previous reviews. In particular, we will show that the justification 
of the procedure is in fact insensitive to the appearance of "resonances" which can arise in 
the multi-phase case, which is the basis for its widespread use in the multiphase situation. In 
all cases we will consider the regular Whitham system defined above without any assumptions 
about the form of higher corrections to the main approximation (1.5) of a slowly modulated 
solution. 

In this chapter we consider general constructions connected with the Dubrovin-Novikov pro- 
cedure and prove some technical lemmas needed for its justification. At the end of the chapter, 
we give a detailed justification of the procedure in the single-phase case. In the next chapter we 
will consider multiphase case, where the justification of the Dubrovin-Novikov procedure will 
have a more complicated form. In conclusion, we consider a simple but typical example of the 
justification of the Dubrovin-Novikov procedure using the Gardner-Zakharov-Faddeev bracket 
for KdV. 

As we said above, we will now consider the general structures connected with the construction 
of the Hamiltonian structure for the regular Whitham system. 

Let us note that according to (2.19) the flows 

^\.=S^^{^,^,,...) = W\x),n (2.23) 

generated by the functionals I'^ according to the bracket (2.16) commute with the initial 
flow (2.3). The flows (2.23) leave invariant the full families of m-phase solutions of (2.3) as 
well as the values U^ = I^ of the functionals I'^ on them. For a complete regular family 
of ?7i-phase solutions with independent parameters {k^, . . . , A:™') it's not difficult to show that 
the flows (2.23) generate linear (in time) shifts of the phases Oq with some constant frequen- 
cies 00°"^ {\J), such that 

5*'^(*,A;^*e^,...) =6<j"'^(U)$^„(6>,U). (2.24) 



According to Lemma 2.1 we have also that the functionals A;" (I) should generate the zero 
flows on the corresponding family of m-phase solutions of (2.3). For V^ coinciding with the 
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values of /^ on A we get then the relations 

'(U) = 0, a,(3 = l,...,m. (2.25) 






Let us give here the definition of a regular Hamiltonian family of m-phase solutions of sys- 
tem (2.3) and of a complete Hamiltonian set of the functionals (2.18). 

Definition 2.1. We call family A of ?7i-phase solutions of system (2.3) a regular Hamiltonian 
family if: 

1) it represents a complete regular family of m-phase solutions of (2.3) in the sense of 
Definition 1.1; 

2) the corresponding bracket (2.16) has on A constant number of annihilators A'^"^, . . . , N'^ with 
linearly independent variational derivatives 6N /6ip^{x) which coincides with the number 
of independent annihilators in the neighborhood of A. 

Let us say that according to the generalized Darboux theorem we can identify the num- 
ber of the variational derivatives 5N /5if^{x) on A with the number of linearly independent 
quasiperiodic solutions v^ (x) of the equation 



Y,Bll^{'P,<Px,..-) 



'{k) 

k>0 



v^'^ -0 

'^j,kx — "-*' 



where v^ (x) have the same quasiperiods as the corresponding functions ip{x) on A. 

Definition 2.2. We call a set (/^, . . . , / ) of commuting functionals (2.18) a complete Hamil- 
tonian set on a regular Hamiltonian family A of ?TT,-phase solutions of system (2.3) if: 

1) the restriction of the functionals (/^, . . . ,/ ) on the quasiperiodic solutions of the family 
A gives a complete set of parameters {U^, ■ ■ ■ , U^) on this family; 

2) the Hamiltonian flows generated by (/^, . . . , / ) generate on A linear phase shifts of Oq 
with frequencies a;^(U), such that 

rank||a;°''(U)|| = m; 

3) the linear space generated by the variational derivatives 61'^ /5ip^{x) on A contains the 
variational derivatives of all the annihilators N'' of the bracket (2.16), such that 



6N^ 



5ip'^{x) 



^7i(U) "' 






5c^*(x) 



for some smooth functions 7f/(U) on the family A. 

Let us note that it follows from Definition 2.2 that if a complete Hamiltonian set of integrals 
(I^, . . . , / ) exists for a regular Hamiltonian family A then the number of the additional para- 
meters (n^, . . . , n*) discussed above is equal to the number of annihilators of the bracket (2.16). 
Indeed, according to Definitions 2.1 and 2.2, the number of the functionals I^ having linearly 
independent variational derivatives on A exactly equals to m + s, where s is the number of 
annihilators of the bracket (2.16). The total number of independent parameters U^ on A is then 
equal to 2m + s due to the wave vectors A:", a = 1, . . . ,m, which implies the above assertion. 

As follows from the condition (2) of Definition 2.2 and from the invariance of the functio- 
nals A^ and I"^ with respect to the flows (2.23), the values 7f,(U) can always be chosen indepen- 
dent on the initial phase shifts on the family A. The values 5/'^/5c^*(x)|a are linearly dependent 
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on A, so it's natural to choose a complete linearly independent subsystem. Remembering that 
the variational derivatives of J'^ (2-7) are linear combinations of the regular left eigenvectors 
jfui^^ + 6q), we can write 



5N^ 



S(p^{x) 



= Y,4{U)K^ljAk{\J)x + eo) (2.26) 

<p(x)=*(k(U)x+eo,U) q 



for some smooth functions nUXJ). The functions ?T-i(U) are then uniquely determined on A and 
we have rank||n' (U)|| = s by Definition 2.1. 

We will need now the construction of the Dirac restriction of a Poisson bracket on a subma- 
nifold. We describe here briefly this construction. Using the terminology of finite-dimensional 
spaces we can say that the Dirac restriction of a Poisson bracket on a submanifold J\f C Ai"' is 
associated with a special choice of coordinates in the neighborhood of the submanifold A^^. The 
coordinates in the neighborhood of the submanifold M are divided into the "coordinates on the 
submanifold" (C/^, . . . , U^) and the constraints {g^, . . . , g^~^) which define the submanifold M^ . 
It is assumed that the submanifold M is defined by the conditions 

5*(x) = 0, i = l,...,n- /c, 

while the functions C/^(x), . . . , [/ (x) on M^'' play the role of a coordinate system on J\f after 
the restriction to this submanifold. 

If the Hamiltonian flows generated by the functions U^{'x) leave the submanifold A^ invariant, 
i.e. we have 

{C/^(x),5'(x)} = for g(x)=0, 

then the pairwise Poisson brackets of the functions U^x.) define the Poisson tensor in the 
coordinates {U^ , . . . ,U^) after restriction on J\f^ which is called the Dirac restriction of the 
Poisson bracket {•, •} on the submanifold J\f G M.^. 

In general, according to the procedure of Dirac, if we have some constraints 5*(x) which 
define the submanifold J\f and some functions U^ (x) which give a coordinate system on J\f we 
must find k linear combinations /3s(U)5'^(x) at every point of AA^, such that for the functions 

U^x) = Wi^) + /3f (U)/(x), j = l,...,k, 

we have the relations 

{c7^(x),5Xx)}=0 

for g(x) = 0. 

The functions f/-'(x) take the same values that U^(x.) at the points of M^ and we can define 
the Dirac bracket {■,-}d on J\f by the formula 

The functions /3s (U) are determined from the linear systems 

{5^(x),5^(x)}|_^,/3^(U) + {5'(x),t/^(x)}|_^, =0, i = l,...,n-k, 

and we can also write 

{U\W}^ = {C/^(x),[/^(x)}|^, -/3l(U){9^(x),5''(x)}|^,/3^-(U) 

for the Dirac bracket on Af^ . 
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We now describe the procedure of the Dirac restriction, which wih be needed in our situation. 
We consider now system (2.3) which is Hamiltonian with respect to some local bracket (2.16) 
with a local Hamiltonian function of the form (2.17). We first introduce the extended space of 
fields 

<^{x) ^(f{e,x), 

where the functions (p{0,x) are 27r-periodic with respect to each O"', and define the extended 
Poisson bracket 

k>0 

We also make the replacement x — >• X = ex and define the Poisson bracket 

{^\e,X),^^{e',Y)}=^e'B'^^^{^,e^x,...)6^'HX-Y)5{0-e') (2.27) 

fc>0 

on the space of fields ip{0, X). 

Let us define again the submanifold fC by the formula (2.12) in the extended space of fields. 
We thus assume that the functions <p{0, X) represent functions from the family A of the m-phase 
solutions of (2.3) with some parameters U(X) at every X. 

It is convenient to choose the boundary conditions in the form U(X) — t- Uq, X — )• ±00, such 
that A;"(Uo) = and define the functions S^^X) by the formula 

r"+00 



1 r+co 

5"(x) = - / sgn(x - y)A;"(y)dy. 

2 J~oo 



We can write then 

ip\9,X) = ¥( ^^'^^^^^ +6>,U(X)), (2.28) 



where ^"[U] is a single- valued functional of coordinates U(X) on the submanifold /C.^ 
Let us introduce the functionals J'^{X) on the functions (p{6,X) by the formula 



J''(X) = l ■■■I P''i^^^'fx,e'^xx,...)j^:^ (2.29) 



and consider their values on the functions of the family /C. 
We can write on /C, 



r{X) = U^{X) + Y.''JiiM)' (2.30) 



1>1 



where JY^JX) are polynomials in the derivatives Ux, UxXj • • • with coefficients depending on U 
and have grading degree I in terms of total number of derivations with respect to X. 

The transformation (2.30) can be inverted as a formal series in e, such that we can write 

U''iX) = r{X) + Y,e'U^i)iX) (2.31) 



^The case when the values k" = are absent on the space of the parameters U requires just a simple 
modification of the definition of S{X) which we do not consider here. 
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on the functions of the submanifold fC. In the formula (2.31) the functions UYj^. are functions of 
J, Jx, Jxx, • • • 1 polynomial in the derivatives Jjc, Jxx, • • • , and having degree I in terms of the 
number of derivations w.r.t. X. 

The values of the functionals J^{X) or U^{X) on the functions from fC can thus be chosen as 
coordinates on /C. We can define also the functionals U(X) on the whole functional space using 
the definition of the functionals J(^) and relations (2.31). Let us introduce also the constraints 
g'^{9^ X) defining the submanifold /C by the conditions g^{9, X) = and numbered by the values 
of e and X 

g\e, X) = ^\e, X) - ¥ /"^MOMl + e, v[3]{x)] . (2.32) 

The constraints g^{0,X) are functionals on the whole extended space of fields (p'^{9,X) by 
virtue of the corresponding definition of the functionals J'^{X). 

The constraints (2.32) are not independent since the following relations hold identically for 
the "gradients" 6g'{0,X)/6ip^{e',Y) on the submanifold /C 



oo /■27r r2Tr 



oo Jo 



5r{Z) 



V(^,^) 



8g\e,X) 



^6^^{e',Y) 



^p^lX.O. (2.33) 



Nevertheless, it will be convenient here not to choose an independent subsystem from (2.32) 
and keep the system of constraints in the form (2.32) keeping in mind the existence of identi- 
ties (2.33). 

Thus, we consider now the values of the functionals [U(X), g(0,X)] as a coordinate system 
in the neighborhood of tC with the relations (2.33). The values of the functionals U(X) will be 
considered as a coordinate system on /C. 

Using the explicit expression for the quantities 5J^{Z)/5ip^{6^X) on /C 



5r{Z) 



K 



5ip^{e,x) 

= Enr(o(*(^ + ^>U(Z)),eA$(^ + 0,U(Z)),...)e'5W(Z-X) 
we can write the convolution (2.33) in the form of action of the operator 



rf^Enr,„(K^-.H'4K^-.H'->'^<-' 



on the distributions 5g'{e,Z)/6if>>{e' ,Y)\k:. 

Pairwise Poisson brackets of the functionals J'^(X), J^^{Y) have the form 



••• / Al^{^{e,X),e^x{e,Xl...)——eH^'\X-Y), 

where 

^7((^, eipx, . . . ) = ^dxQ'''\<^, eifx, ■■■)■ (2.35) 

Substituting the functions if^{6,X) in the form (2.28) it is easy to see that the Dubrovin- 
Novikov bracket is the main part (in e) of the pairwise Poisson brackets {J'^{X), J^{Y)} on the 
submanifold /C in the coordinates [U(Z)]. 
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The Poisson brackets of the fields ip^{6,X) with the functionals J^{Y) can be written as 
W\e,X),J^{Y)} = Y,^'Cl^iv>{0,X),e^x{e,X),...)6'~^HX-Y) (2.36) 

l>0 

with some smooth functions Cl^Jif, €(px, ■ ■ ■)■ 
We also have in this case 

Cl!^^iipie,X),e^x{0,X),...)^S'^{^{e,X),e^x{0,X),...) (2.37) 

by virtue of (2.23). 

For any function of the slow variable q(Y) we can write 

\ip\e,X), f qiY)J'^{Y)dY\ =Y,^'Cl';^i^iO,X),e^x{0,X),...)qix. (2.38) 

The leading term in the expression (2.38) on fC has the form 



^\e,X), j qiY)J^iY)dY 



KH 



^(0) (* (^ + ^' ^) ' k^^e^ (^ + e,xy...] q{X) 

q{x)s^^ U{^ + e,xY k^^,, (^ + o,x 



5 • • • I 5 



where S^^{(f, cpx, • • • ) is the flow (2.23) generated by the functional I^. According to (2.24) we 
can write then 



ipHe,X),Jq{Y)J^iY)dY 



a;°'^(X)$^. ( ^^ + 6, U(X) ) q{X) + 0(e). (2.39) 



/c 



Similarly, for any smooth function Q{6, X), 27r-periodic in each 9°', we can write on the basis 
of (2.36) 



jj ... jj Q(^^ + e,x^ y\e,x),j>^{Y)] 



\1C (27r)™ 
r.27r r2TT 



/•2-iT 1-Ztt rl"^/? 

= u'^^{X) ... Q(0,X)$^40,U(X))— — <^(X-y)+0(6). (2.40) 

JO JO l^^j 

In view of relations (2.30), (2.31) for the functionals J^^{Y), U^{Y) we can also write 
If'iO.X), j q{Y)U"{Y)AY 

...j^ Q{e,x)^^{eMx))j^-y^5{x-Y) + o{e). (2.42) 

In this case, by virtue of (2.25) we have 



uj^^iXWe. (^ ^ + e, U(X) ) q{X) + 0{e), (2.41) 



ip\e,X),JqiY)k''{V{Y))dY 



= 0(e), (2.43) 

fC 
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i-2tt i-2-n / ^( Y\ \ (\"^f) 

I ■■■I Q(^^ + 0,XJ{,,^(0,X),^"(U(y))}|^^^^ = O(6) (2.44) 

for fixed values of coordinates [U(Z)]. 

Let us note here that since we consider the change of coordinates, which depends exphcitly 
on e, we shah assume now, that aU relations are written in the coordinate system [U(X), g(0, X)] 
near /C. In particular, all the values on /C will be functionals of U(X), represented usually in 
the form of graded expansions in e with coefficients depending on U, Ux, Uxx, • • • • 

Let us prove here some lemmas about the structure of Poisson brackets on the submanifold /C 
which we will need in the further consideration. 

Lemma 2.3. Let the values U^ of the functionals I^ on a complete regular family A of m-phase 
solutions of (2.3) be functionally independent and give a complete set of parameters on A, ex- 
cluding the initial phases. Then for the Poisson brackets of the functionals k"(U{X)) and J^{Y) 
on /C we have the following relation 

{k^{U{X)),J^{Y)}\^ = e[u:"^{\J{X))6{X-Y)]^ + 0{e'). (2.45) 

Proof. The conditions of Lemma 2.3 coincide with the conditions of Lemmas 2.1 and 2.2. 
Consider the Hamiltonian flow generated by the functional J q{Y)J^^(Y)dY according to the 
bracket (2.27) with a compactly supported function q{Y) of the slow variable Y = ey. 
According to (2.38) we can write at the "points" of the submanifold /C 



"Pt 



q{X)u:P^{\^{X))%, {^+eMX)) + er/j^, (^ + 0,x) +0{e^) (2.46) 



with some (27r-periodic in each 0") functions rf, A6,X). 

Let us introduce the functionals fe"(J(X)) = k°^{J^{X), . . . ,J^ {X)) using the same func- 
tions /c°(U). According to (2.30), (2.31) we can write 

{k^{\J{X)),J^{Y)]\^ = {k^{3{X)),J^{Y)]\^ + 0{^). 

Consider the evolution of the functionals A;"(J(X)) according to the flow generated by 
J(7(y) J''(y)dy on the submanifold /C. Using relations (2.46) we can write 

r)k'^ f)k" 

kf{J{X)) = ^^ (J(X)) jr(X) = ^^ (J(X)) q{X)oo^'^{V{X)) 



f...fEnr.„(K^--<^-')'4K^--(-' 



J d' ^r ^S(X) , ,,,,^A d™^ 



^^'i^%A^T^ + (^Mx) 



dXi ^'^ V e ' ^-^ 'J (2vr)^ 
— (U(X)) {qiX)u^^iViX)))^ 



r2iT (■2-K 



\^ ■■■l^Y. ^^m (* (^' u(^)) , ^'^r' {0, u(x)) , . . . ) 

Jo Jo ;>i 

X fc^H^) • • • k^'--{X)<^l,,^,,„,.,_, {0,V{X)) 



(27r)™ 

aha /•2tt f-2n 

'u(x)) / • • • / E nrw (* (^' u(x)) , F*e. {e, v{x)) ,... 
Jo Jo r^ " 



^^dU 

l>0 



(F'e 



xk-^iX)...k'^^iX)r^l^^,^,,„,.,ie,X)j^^ + 0{e'). 
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It's not difficult to see that the first part of the above expression contains the integrals over 6 
of the expressions P!^0 {6, X) and is equal to zero. It is easy to see also after integration by parts 
that the third part of the above expression represents the value 



e 



2n f2TT OTQ ^ma 

.0 aZ7^(U(^»4u(x),W4](^,^)(^ 



and is equal to zero according to Lemma 2.1. 

Thus, we can write in the main order the expressions for the evolution of the functionals 
k'^{J{X)) on the submanifold /C in the form 

knJ{X))=e{q{X)u^^{\JiX)))^ 

piUa /•2n f2n Amo 

Using relations (2.13) we get then 

fcf (J(X)) = e [q{X)co"^{V{X))]^ + 0{e'), 
i.e. 

{k'^iJiX)), J''{Y)}\jc = eu'''^{V{X))6'{X -Y) + euj''/6{X -Y) + 0{e'), 
which implies (2.45). ■ 

Easy to see that according to Lemma 2.3 and relations (2.30), (2.31) we can write also 
{A;"(U(X)), U'^{Y)}\k = e [u'^^{V{X))6{X - Y)]^ + 0{e^) (2.47) 

on the submanifold /C. 
Corollary 2.1. Under the conditions of Lemma 2.3, we have also 

{k-{X),k^iY)}\^ = 0{e') (2.48) 

for the functionals k(U(X)). 

Indeed, by virtue of (2.25) we have 



{k-{X),k^iY)}\^ = {k-{X),U^{Y)}\^ — {U{Y)) 

= eu^'^{X)kUx)6'iX -Y) + ,oj'^^{X){kUx))j{X - Y) 
+ 6 (^"^(X))^ kUx)5{X -Y) + 0{e') = 0{e'). 

Let us note now that in the proof of Lemma 2.3, we have not used the fact that the function- 
als J'^{Y) belong to our special set of functionals (2.29) and used only the fact that the flow 
generated by the functional I^ leaves invariant the family A generating linear shifts of 9q with 
constant frequencies a;"^(U). We can therefore formulate here the following lemma. 

Lemma 2.3'. Let the values V^ of the functionals I'^ on a complete regular family A of m-phase 
solutions of system (2.3) be functionally independent and give a complete set of parameters on A, 
excluding the initial phases. Let the flow generated by the functional 



1= I P{^,^„...)dx (2.49) 

leave invariant the family A generating a linear shift of 6q with constant frequencies a;°(U). 
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Consider the junctionals 

r2n f2Tr 



J{X) 



P(ifi,e(^x,(^ i^xx, 



dre 



(2.50) 



/o JO '(2^)™' 

Then for the Poisson brackets of the Junctionals k°'{\]{X)) and J{Y) on fC we have the relation 
{k^mX)l J{Y)]\^ = e [(I;"(U(X))<5(X - Y)]^ + 0{e^). 
Proof of Lemma 2.3' completely repeats the proof of Lemma 2.3. 

Lemma 2.4. Let the values V^ of the functionals I'^ on a complete regular family A of m-phase 
solutions of system (2.3) he functionally independent and give a complete set of parameters on A, 
excluding the initial phases. Then for the constraints g^{6,X) imposed by (2.32) and smooth 
compactly supported function q{X) as well as smooth 2-K-periodic in each 6°' function Q{0,X) 
we have the following relations on the submanifold K, 



g\e,X), j q{Y)J>^{Y)AY^ 

S{X) 



0{e) 



Q 



+ e,X){g\0,X),J^^{Y)} 



(T^e 



UO JO \ - J (2vr)^ 

Proof. Indeed, by (2.39), (2.40), and Lemma 2.3 we have 

'S(X) 



0{e). 



(2.51) 



K. 



g\e,X), q{Y)J^^{Y)dY 



^l^ 



/Cfi 



+ e,\J{X)\uj'"'{X)q{X) 



- l^^ (^ + ^' U(^) ) / sgn(X - Z) {u:^^{Z)q{Z))z dZ = 0, 



27r 



UO 







'\i^^ + e,x]{g\e,x),j^{Y)] ^"^ 



{2^y 



uj''^{X) 

r2Tv 



2-K 



e 

2-K rl"^fl 

Q{e, X)%^ {6, \^{X))——5{X - Y) 



m 







JO (2vr)- 

2iT r\^^R 1 

Q(0,X)<l>^.(0,U(X))^,^-^-sgn(X - Z) {uj'-^{^Z)b{Z - Y)\^dZ ^ 0. 







(27r)'" 2 



Similarly to the previous case, we can formulate also the following lemma. 

Lemma 2.4'. Let the values U^ of the functionals l^ on a complete regular family A of m-phase 
solutions of system (2.3) be functionally independent and give a complete set of parameters on A, 
excluding the initial phases. Let the flow generated by the functional (2.49) leave invariant 
the family A generating a linear shift of 9q with constant frequencies a}"(U). Then for the 
constraints g^ {6 , X) and the functionals J{X) imposed by (2.50) we have the following relations 
on the submanifold fC 



g\e,X), q{Y).J{Y)dY 



Oie), 



K. 
2^ r2^ /g(^) „ „\ r ,-,„ ,.. ;._.. d'"^ 



Q 



+ e,x\{g\e,x),j{Y)] 



(27,y 



0(6). 



K. 



Let us consider now the Dirac restriction of the bracket (2.27) on the submanifold /C. 
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For the Dirac restriction of the bracket (2.27) on the submanifold /C we have to modify the 
functionals U^{Y) by Hnear combinations of the constraints g^{6,X) 

„,'S[U](X) , ^ ,^^^ ^^^\ d™0 



U^'iY) -- 


J-oo Jo 


Jo 



g\e,x)a'l ^-^^ + e,[v],x,Y j j^-y^dx, 

so that the flows generated by the functionals U^(Y) leave /C invariant and then use the func- 
tionals U^(Y) to construct the Dirac bracket on /C. 

Since the functionals U^{Y) are defined with the aid of the functionals J{Z), technically it 
is more convenient to modify the functionals J^(Y) 

J>(y) = J^iY) + p_^j^^. . . j\\e,X) ;3m|^ S[U[J]](X) ^ ^ ^^^^^^^^^ ^^j ^dX, 
so that the flows generated by the functionals J^{Y) leave /C invariant and then put 

?7^(y) = j>(y) + j;e't/5)(y) 

just using the functionals iiY) instead of i{Y) in (2.31). Both the approaches give, certainly, 
the same result for the Dirac restriction of the bracket (2.27) on /C. 
The distributions /3f (S(X)/e + 6,X,Y) must satisfy the relation 

/+00 j-2-K j-2tt / ^( 7\ \ r\"'-ff' 

.J.-L W<'.A-).^«''.z)}U?(4)..'.z.r)^dz 

+ {g\e,X),J^{Y)]\^ = Q (2.52) 

on fC and are defined at each "point" of /C up to linear combinations of the distributions 
8J^{W)/5ip^{6,X) by virtue of the original dependence (2.33) of the constraints (2.32). 

To obtain a local Poisson bracket on /C we are trying to find the functions /3f (0, X, Y) in the 
form 

/3f (0, X,Y) = Y, e'/3f(;) [e, X, Y), (2.53) 



where the functions (3'^,j^Jd,X,Y) represent the local distributions 

p=0 

with the total grading / in e, considering that the derivatives of the delta function 6^^'{X — Y) 
have by definition the degree p. 

Thus, we assume that all the functions /3^yz) »(^' -^) ^^ ^ ^^^ local functionals of (U(X), Ux, 
. . . ) at every 0, polynomial in the derivatives (Ux, • • • ) and having degree / — p under the 
previous definition. Such a structure of I3^{6,X,Y) is obviously equivalent to saying that the 
functionals 

/ + 00 
J^^iY)q,{Y)dY 
-oo 

can be modified by linear combinations of constraints (2.32) with the coefficients 

l>l l>l p=0 
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so that the corresponding flows leave invariant the submanifold /C. Under this scheme, the 
derivatives d^q^/dX^ of slowly varying functions q^{X) also have grading p as the derivatives 
of the functionals U(X). 

The pairwise Poisson brackets of the constraints g^{6,X), g^{0\Y) on /C can be written in 
the form 

{g\G,X),g\e\Y)]\^ 

= W{e,x),^{e\Y)]\K,- W{e,x),u\Y)]\K.%,{^ + e' ,\]{Y) 

-^\j^ {^^ + eMx)^{u''{x),^{e\Y)]\^ 

+ cD^j,. (^ + 6, u(x)) {u-{x), u\Y)}\^ %, {^ + e\ u(y) 

- ^ jw{e, X), k^{w)}\,c sgn(y - w)dw '^l„ (^ + e', u(y) 

- Y^^lc. (^ + 0, u(x)) I sgn(x - z){fc-(z), (^^(0', y)}kdz 



+ ^^^c(^ + 0,U(X) 



x|sgn(X-Z){A;-(Z),^^(y)}kdZcl>{^, (^ + 0',U(y)) 

X l{U''{X),k''{W)}\,csgn{Y-W)dW^i,J^ + e',V{Y) 
+ ^n.(^^ + e,V{X)^ (2.54) 

X f sgn{X - Z){k''{Z),k''{W)}\!csgn{Y -W)dZdW'^i„ (^^ + 0' ,\J{Y)] . 

For the purposes of this chapter, we need only "approximate" restriction of the bracket (2.27) 
on the submanifold /C. Specifically, we will need to prove the existence of only the first term 
(3'^,-^Jd,X,Y) in the expansion (2.53) having the form 

/3f(^)(0, X, Y) = /3f(^)_^(0, \J{X))5'{X -Y)+ /3f(,)^o,(0, \J{X))U^5{X - Y), (2.55) 

which is equivalent to the approximation 

A[q](i)(0,X) = /3f(,) ,(0,U(X))g^x(X) + /3f(^)o,(0,U(X))C/^<ZM(^) (2-56) 

for the functions Bi[q^{6, X). 

To find it we need only the leading term of the convolution of the Poisson brackets of con- 
straints on /C with i?jrqiM)(S(y)/e+0',y) and the first (in e) non- vanishing term of the brackets 
of constraints and the functional Jrqi . 

Note that according to the relations (2.43), (2.44), (2.47), and (2.48) all the terms of the 
bracket (2.54) behave as the values of the order of at most 0(1) for e — )■ in the convolution 

withi?,[q](i)(s(y)/e + 0',y). 

We impose additional conditions 

J^ ■■■j^ <^U^,ViX))P^^^^ie,X,Y)j^^O, a = l,...,m, (2.57) 
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or equivalently 

r-27r /.27r 



/ •••/ $^,40,U(X))S,[q](i)(0,X) 

Jo Jo 



dre 

(27r)™ 



0, a = 1, , 



(2.58) 



to be confirmed aposteriori for the functions /3^.^s(0,X, y) and Bj[q^(^i^{6, X). The Poisson 
brackets of constraints (2.54) on /C can then be reduced to the form 



{g\0, X),g^ {e\ y)}k ^ {g\e, X),g^ {6', Y)}\f 

= W\e,x),ip^{e',Y)}\,c-W\0,x),uHY)}\K^l 



u^ 



s(y) 



+ 6>',u(y) 



^U" 



1 



e 

S(X) 



+ eMx)]{u''{x),^\e',Y)}\,c 



2e 
1 



^^c 



+ ^^^^ 



e 

e 

S(X) 



+ e,vix)]{u''{x),u\Y)}\K^' 






s(y) 



+ 0',u(y) 



+ e, u(x)J y sgn(x - z){A;"(z), ^^(0', y)}kdz 

+ 0, U(X)) / sgn(X - Z){k^{Z), U^{Y)}\jcdZ 



^Ui^^ + ^'^'^O')] 






(2.59) 



in finding 5j[q](i)(0,X). 

To find the leading term of the convolution of the Poisson brackets of constraints on /C and 
i3jjq](i)(S(y)/e + 0',y) we can simplify further the relation (2.59). Namely, taking into account 
relations (2.41), (2.47) we can omit in the leading order the second and the last terms of (2.59) 
in the convolution with B jy^^i-^{S{Y) / e + 0' , Y). By (2.42) and (2.58), we can omit in the leading 
order also the third term of (2.59) in the convolution with B jl^ii\{S(Y) / e + 0' ,Y). The fourth 
term of (2.59), obviously, is of the order 0(e) and also can be omitted in the leading order. 

As a result, the principal term of the convolution of (2.54) with the functions 
e5j[q](]^)(S(y)/e + 0',y) can be written in the form 



{5^(0, X), 5^(0', y)}kei?,[q](i) [^-^ + e\Y\ .:^— dy 



X fc"H^)---fc"^(^)i?,[q](l)e'^l...e". (^ 

^+0,U(X))|sgn(X-Z) 



(2^)' 
S{X) 



[1] 
+ 6>,U(X) 



+ e,x 



2 '' 



Jmz),^io'y)}l..B,,,,m.e,Y 



(27r)"^ 



-dy 



dZ. 



[2] 



(2.60) 



Let us remind here that the indexes [1], [2] mean as before the terms of the corresponding 
graded expansion on /C having degree 1 and 2 respectively. 
We have also 



{g\e,x),j[^]}\jcii} = W\e,x),j[^]}\jcii]-n 



S{X) 



+ e,ij{x) {u^x),j^^]}\icii] 
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+ e,V{X)^ J sgn{X - Z){k"{Z),J[^]}\lC[2]dZ. 



(2.61) 



Let us prove here the following lemma which allows to investigate the further construction of 
the functions i?j[q](i)(0,X). 

Lemma 2.5. Let the functions Biiq]fi\{0,X) satisfy the conditions (2.58), i.e. 
and the relation 







M(X)S,[q](i)(0,X) + 4]fq](0,X) = 0, 



where 



BUX) 



(2.62) 
(2.63) 



QS 



Y,Bl^meMx)),k%x)^e'^{e,\j{x)),...)k^-{x)...k'^^{x)-—-_ 



A 



[i][q] 



s>0 

S(X) 



de^'s ' 



+ e,x 



(2.64) 



{V9*(0,X),J[q]}|^[l]-C|.[;. 



S(X) 



+ 0,U(X) {[/'^(X),J[q]}|^[l]. 



T/ien the functions i?jrqi(i)(0,X) satisfy the relation 



{g^(0,X),5^(0',Z)}k6S,[q](l) ( ^+0',z) ^^^ 



[1] 



Proof. The fulfillment of condition (2.58) allows to reduce the expression 

{5^(0, X), 5^(0', z)}k6i?,[q](i) (^^ + e',z\ S^dz 



(2.65) 



e J (27r)''" 

according to formula (2.60) up to the terms 0{e^). From (2.62) we can get then 

{g^{e,X),g^e',Z)}\^eB^^^]^,^(^ + e',z) dZ + {g\0,X),J^^]}\,c 



{27Ty 



2^^n e 



+ e,\JiX)] Jsgn{X-Z) 



{A;-(z),,.^(0My)}kei?,[q](i) (^ + e',w^ ^dW 



dZ 



[2] 



W,^ ( ^^ + e, V{X) ) I sgn(X - Z){k"iZ), J[q]}k[2]dZ + 0(e2). 



(2.66) 



The proof of the Lemma follows then from the orthogonality conditions of the left-hand part 
of (2.66) to the variational derivatives 



Skl^{3{X')) 



6ifie,x) 



K 



dk\ , 5J\X') 



13 = 1, . . . ,m 



(2.67) 



JC 
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according to (2.33) in the order e^ on /C. Here we use the functions k^(\J) to define the 
functionals k"{3{X')). Let us remind also that the functionals U(X) and J{X) are connected 
by the transformations (2.30) and (2.31). 

Indeed, the right-hand part of (2.66) has the form 



n. (^ + e,lJ{X)] ysgn(X-Z)Qfqj(Z)dZ + 0(e2) (2.68) 



with some functions Q?^{Z). 

The convolution with the quantities (2.67) can be represented as the action of the operator 

on the functions of S(X')/e + 9 and X' . 

Expanding the action of the operator (2.69) in powers of e we can write it's main part in the 
form 

„ ■■■/ (ifixW^')) 

x^n>„(^(...),fi.».(...)....)fc'"(xO---t-(x-) ^„,^^' 

l>0 

We can see then after the integration by parts that the main part of the operator (2.69) 
reduces to the convolution (w.r.t. 6) with the vectors 



0^ (U(^ )) C,[u(x')] [-^ + ^) 



which are identically equal to zero according to Lemma 2.1. 

We can conclude, therefore, that terms of the order of O(e^) in the right-hand part of (2.66) 
do not play any role in the convolution with the values (2.67) in the order e^. 

Let us consider now the action of the operator (2.69) on the terms of the order of e in the 
expression (2.68). We can see that the terms in which the operators d^/dX"* apply only to the 
functions *0c(S(X')/e + 6, U(X')) yield the quantity 






identically equal to zero. 

The terms of the order of e^ occur at the same time with a single differentiation of the function 
sgn(X' — Z) and application of the remaining derivations d/dX' to ^ga (S(X')/e + 0, \]{X')) in 
the leading order. In the required approximation the functionals i{X) should be also replaced 
by U(X) in the argument of dk^ /dU'^. 

The corresponding values are then equal to 






r,^ (U(X')) 

f2TT r2lT 
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Using again the relations (2.13), we obtain that the orthogonality of the right-hand part 
of (2.66) to the values (2.67) in the order e^ is equivalent to the relations Q^ AX') = 0. Given 
that the main part of (2.68) represents the difference between the left-hand sides of (2.65) 
and (2.62) we obtain the assertion of the lemma. ■ 

For a regular Hamiltonian family A and a complete Hamiltonian set of integrals (I^, . . . , / ) 
we can also prove the following lemma. 

Lemma 2.6. Let the functions Bj\^ii\{0,X) satisfy conditions (2.62). Then the functions 
i3jjq](i)(0,X) automatically satisfy (2.58). 

Proof. We first prove the following statement. 

The values Qmix)]^^) ^^^ orthogonal (for any X) to the values Al^,, A6,X), i.e. 

2tt i'2-k c\"^B 

CTU.W4ilfal(^>^)7T^Z^ = 0. (2.70) 



Indeed, as we know, the values {(7*(0,X), J[q]}|A: identically vanish under action of the op- 
erator (2.34) (with the replacement of Z to X). In the leading order (~ e) it's action is given 
by the convolution w.r.t. Q of the values {g^{0, X), J[q]}|A:[i] , imposed by (2.61), with the values 
Ql^jj,-^.JS{X)/e + 6). We know also that the values Q^-^,-^^AS{X)/e + 6) are automatically 

orthogonal to the functions $ga(S(X)/e -|- 0,U(X)), so we get the relation (2.70). 
Thus, the implementation of (2.62) implies the conditions 



/'^•••/^iu(x)]W^g]W^.M(i)(^'^)- 



d'"^ 



(27r)'" 
which is equivalent to 

j'2-K /•2tt c\^^f) 

w"'^(U(X))y^ ■••j^ $^^40,U(X))i?,[q](i)(0,X)^^-^ = O 

view the skew-symmetry of BlL{X). 

From condition (2) of the definition of a complete Hamiltonian family of commuting func- 
tionals we immediately obtain now conditions (2.58). ■ 

It is not difficult to see that from the relation (2.70) follow also the conditions 
r2-K f2T[ (V^ft 

'^Ix-.m^m.^i^^X)—-^^ (2.71) 



/•ZTT pi 

Jo Jo 



\[lJ(X)]'^'')^llM^''^^)J^r 



by Proposition 2.1. 

We come, therefore, to investigation of system (2.62) to construct an "approximate" restric- 
tion of the Poisson bracket (2.27) on the submanifold fC. In the remaining part of the article the 
investigation of system (2.62) will play the basic role for the justification of the main results. 

In what follows we consider separately the single-phase (m = 1) and the multiphase (m > 2) 
cases. 

The following lemma can be formulated for the single-phase case m = 1. 

Lemma 2.7. Let A be a regular Hamiltonian family of single-phase solutions of (2.3) and 
(I^, . . . ,/ ) be a complete Hamiltonian set of the first integrals of the form (2.18). Then the 
functions Bji„ui-^{d,X) can be found from system (2.62) and can be written in the form (2.56) 
with smooth dependence on parameters \J{X). 
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Proof. System (2.62) in the single-phase case is a system of ordinary differential equations 
in 9 with skew-symmetric operator nlAX). It is easy to see also that the right-hand side of 
system (2.62) has the form 

r {e, U(X)) q,,x{X) + vT (0, U(X)) U^q,{X) 

with periodic in 9 functions e^(e,U(X)), r]f{9,V{X)). 

The orthogonality conditions (2.71) imply the orthogonality of both the sets of functions 
C>'i9,lJ{X)) and 7]'^{9,\J{X)) to the functions k|['u(X)](^)' ^^^^ ^^^^ system (2.62) can be split 
into independent inhomogeneous systems 

^[0] i^)f^m,i (^' U(X)) = r (0, U(X)) , (2.72) 

^[o] (^)/?Si),OA (^' U(X)) = ry^ {e, U{X)) (2.73) 

defining functions (2.55). 

Both the systems (2.72), (2.73) are systems of ordinary linear differential equations with 
periodic coefficients and the skew-symmetric operator B^,^JX). The zero modes of the opera- 
tor n,^JX) are given by the variational derivatives of annihilators of the bracket (2.16) rep- 
resented as functions of 9 on the manifold of single-phase solutions and are orthogonal to the 
right-hand parts of (2.72) and (2.73) according to (2.26) and (2.71). Eigenfunctions oi n,L{X) 
form a basis in the space of 27r-periodic functions <f{9). Besides that, the nonzero eigenvalues 
of i??Qi {X) are separated from zero in this case. Thus, the 27r-periodic functions /3^^n -^{9, \J{X)), 
I3'^,-^.Q^{9,\J{X)) can be found from systems (2.72), (2.73) up to the variational derivatives of 
the annihilators of the bracket (2.16). If we impose additional conditions of orthogonality of 
l3^f-^-.{9,X,Y) to the variational derivatives of the annihilators of the bracket (2.16) on the 
manifold of single-phase solutions we can suggest a unique procedure of construction of the 
functions l3'^^^^{9,X,Y). The functions /3^^(^)(6',X, y) satisfy (2.52) in the order of 0(e) and 
conditions (2.57) and smoothly depend on the parameters U(X), which implies the required 
properties for the functions i?j[q](i)(0,X). ■ 

Let us formulate now the theorem that gives the justification of the Dubrovin-Novikov pro- 
cedure for a regular Hamiltonian family of single-phase solutions of system (2.3). 

Theorem 2.1. Let A be a regular Hamiltonian family of single-phase solutions of system (2.3). 
Let (J-*^, . . . ,/^) be a complete Hamiltonian set of commuting first integrals of system (2.3) 
on A having the form (2.18). Then the Dubrovin-Novikov procedure gives a Poisson bracket of 
hydrodynamic type on the space of parameters {U'^{X), i^ = 1, . . . , N}, where U^ = {P^) ■ 

Proof. The antisymmetry of the Dubrovin-Novikov bracket follows immediately from the an- 
tisymmetry of the original bracket. Let us prove the existence of the Jacobi identity under the 
conditions of the theorem. 

Consider the Poisson bracket (2.27) and the functionals 

r-27r 



Jl,^{X) = r{X) + ejj^ %^(e,^)/3^i) (^ + e,Z,x) ^dZ 



in the neighborhood of the submanifold /C. (Note that the functions /^^^-i depending on the 
values U(X) on fC, can be considered as functionals in the neighborhood of fC by the definition 
of the functionals U(X)). 

The pairwise brackets of the functionals JY-^AX), J^^-.{Y) can be written as 

{Jl,^{X)J^^^{Y)] = {r{X),J^{Y)] 
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+ el£^ !^riX),g^ie,W)P^^,^(^^+9,W,Y^^^dW (2.74) 

We note here that despite the presence of non-local terms in the brackets of constraints 
{g'{e,Z),g^{e',W)}, and {J''{X),g^ {9,W)}, {g'{e,Z),Jf'{Y)}, the bracket of the functionals 
{Jf^^(X), Jt^-^JY)} on /C has a purely local form because of the conditions (2.57). We can also 
see here that expression (2.74) is a regular at e — t- distribution after the restriction on /C. 

From the form of the constraints and condition (2.57) it is easy to obtain that the pairwise 
Poisson brackets of the functionals Jy^^{X), JttJY) on /C can be written as 

I 
{jfi)(X),J(';)(y)}|^ = j;6'j;aJ^^j(U,Ux,...)^^''^(^->^), (2.75) 

where a'f/^ , are some functions of U, Ux, • • • , polynomial in the derivatives and having deg- 
ree I — q. By virtue of relations (2.51) we can state also that the leading term in e in the 
expression (2.75) is still the same Dubrovin-Novikov bracket on /C, since the corrections given 
by the last three terms of (2.74) are of the order of O(e^) on /C. We thus have 

I 
{j^,^{X)J^^^iY)}\^ = e{U^X),U^{Y)}^^ + Y,e'^a'(^%^{V,lJx,...)S^'HX-Y), 

l>2 q=0 

where 

{U'^iX), c/^(y)}DN = {A7) (U(X)) 6'{X -Y) + -^g^Ul6{X - Y) 
is the Dubrovin-Novikov bracket on the space of fields U(X). Thus, we can write 

+ ^Z74w)^^'^"m(U,Ux,...)'^(^)(X-1^) (2.76) 

^ ^ l>2 q=0 

on the submanifold /C. 

For further discussion it is convenient to introduce "regularized" functionals. Namely, for any 
smooth compactly supported vector-valued function q(X) = {qi{X), . . . ,gjv(X)) we introduce 
the functionals 

J[q] = Jq,iX)riX)dX, J(i)[q] = J q,iX)J^,^iX)dX, U[^^ = j q,{X)U%X)dX. 
We can write then 

hoc i'2-K 



/+00 j-l-K r f SiZ^ \ ■ ^ dO 

+ '£^ r[^^^^^''^^'^^^^l^^W (^+^'^)} 2^d^ (2.77) 
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+ e^ 



Bi[^-^t^,J^ + e,z\g\9,Z),g^{e',W)Bj[p]^,){'^^^^ + e',W] '^^^dZdW. 



( sm 



d9de' 

( 

27r 27r 



We can see again that the last three terms of (2.77) give corrections of the order of e^ on the 
submanifold /C. 

For the expansion of the brackets (2.77) in the neighborhood of the submanifold /C it is 
convenient to use the expressions 



ip'ie, X) = g\e, X) + $^ (^ + e, u(x) 



(2.78) 



making it possible to easily expand any functionals of <^{0, X) with respect to the values of g. 
For example, using the expression (2.20), we can write for the brackets {^[q], <^[p]} in the neigh- 
borhood of K, 



{•^M'-^Ip]} = {'^[q]''^[p]}| 



+ 



2-K 



1 

E 

s>0 



5<fi{e,W) 



J J e'A7 {^{6', X), e^xie', X), . . . ) q,{X)p^MX)^dX 



K 



dO 



where the values { J[q], </[p]}|a: as well as all the other values on /C are calculated on the functions 
(^*(^, W) G /C corresponding to the same values of the functionals [U(Z)] as the original functions 
(/?*(0, VK). We can see that the first term of the expansion of {"^[q]; <^[p]} near /C also has the 
order 0(e) as well as the values {-/[q], <^[p]} on /C view the condition (2.35). 

Considering relations (2.36), we can write a similar representation for the remaining terms 
of the expression (2.77). 

We can see that the last term of the right-hand side of (2.77) has the order O(e^) in the 
whole space which gives the same order in e for the terms of its expansion after the substitution 
of (2.78). 

By Lemma 2.4, the second and the third terms of (2.77) are of the order O(e^) on the 
submanifold /C. However, this can not be said about the next terms of the expansion of these 
brackets in the values of g. To find the values of these terms near /C we write e.g. 



2-K 



B 



i[q](i) 



S{Z) 



dQ 



+ 9,Z\g\0,Z)—dZ 



2tt 



B. 



q](i) 



S{Z) 



d9 



+ e,zW'{o,z)—dz 



2tt jn 

A[q](l)(0,Z)cl>^(0,U(Z))-dZ. 



After rather simple calculations it's not difficult to get then for the second term of (2.77) 
near /C 



" ' A[q](l) (^ + ^' Z) 9'{0, Z), J[p] I ^dZ 



2tt 



II 



" 'Bi^^]^^J^ + 9,Z)g\9,Z),J[p^ 



+ 



B, 



[q]{i) 



e 

S{Z) 



de 
— ( 

27r 



-dZ 



+ e,z 



6Cl^^{^{9,Z),ecpz{e,Z), 



6^p^{e',w) 



xg\e',W)^—dW 



df) 

p^(Z)-dZ + 0{e) 
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+ 



B. 



M]We 



S{Z) 



+ 0,z\^^ sgn(Z - W) {k{W), J[p]} 1^ dW + 0(e) 



X5' 

27r 



Z)^dZ + 0(g2) 







B 



+ 



i[q]{i) 



^i[q](i) 



5(Z) 
e 



+ ^,ZUX^,^),^[F 



k: 



2^' 



dZ 



+ e',z 



dS^nviO',Z),e^ziO',Z),^ 



s^'^{e,w) 



K. 



(\f)' 

p^{Z)—dZ + 0{e) 



>^g\e,w)—dw 



+ 



B. 



»[q]{i) 



S{Z) 



+ e,Z]uj^{Z)p^{Z)+0{e) 



g\e,z)'^dz + o[e). 



In general, considering the expansion of the brackets (2.77) in the neighborhood of /C, we can 
write the relations 



{'^(l)[q]>'^(l)[p]} - {•^(l)[q],«^{l)[p]} 



K. 



+ e 



2-K 



( S{W) 



dQ 



TiHv\ V-^ + ^' W^' A 9\0, W)—dW + 0(g2) 



(2.79) 



with some regular at e — )• 0, 27r-periodic in 9, functions rj[qp](^, W, e). 

Now for the bracket (2.27) we consider at the "points" of /C the Jacobi identity of the form 

{{^(l)[q],^(l)[p]},^(l)[r]}+C.p. = (2.80) 

with some smooth compactly supported q(^), p(^), r(^)- 
We have on /C the relations 

for the bracket of the functionals g^{9, W) and J{i)[r]- 
It's not difficult to get also 

{f/^(W^),J(l)[r]}L = {j^(H^),J(l)[r]}L+0(6^) 

Now using relations (2.76) and expansion (2.79), we can see that the Jacobi identities (2.80) 
on /C coincide in the leading order in e (e^) with the similar Jacobi identities for the Dubrovin- 
Novikov bracket 



'^{f^[q]>f^[p]}DN 



{f/^(iy),[/rrl}DNdT^ + C.p. = 



5W{W) 
on the space of fields U(X). ■ 

Let us formulate here also the Theorems justifying the invariance of the Dubrovin-Novikov 
procedure and the Hamiltonian properties of the Whitham system (2.10) with respect to the 
averaged bracket in the single-phase case. 

Theorem 2.2. Let A be a regular Hamiltonian family of single-phase solutions of (2.3). Let 
{L^, . . . ,1^) and (/'^, . . . , /'^) be two different complete Hamiltonian sets of commuting first 
integrals of (2.3) having the form (2.18). Then the Dubrovin-Novikov brackets obtained using 
the sets {L^, . . . , I^) and {L'^, . . . , I'^) coincide with each other. 
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Theorem 2.3. Let A be a regular Hamiltonian family of single-phase solutions of (2.3). Let 
{I^, . . . ,1 ) be a complete Hamiltonian set of commuting first integrals of (2.3) having the 
form (2.18) and H be the Hamiltonian function for the system (2.3) having the form (2.17). 
Then the Whitham system (2.10) is Hamiltonian with respect to the corresponding Dubrovin- 
Novikov bracket (2.21) with the Hamiltonian function (2.22) 

f+OO 



/+0O 
{Ph) (U(X)) dX 
-oo 



Proofs of Theorems 2.2 and 2.3 are identical for the single-phase and the multiphase cases, 
so we prove them in general m-phase case in the next chapter. 

3 Dubrovin-Novikov bracket and the multiphase case 

We turn now to the general multiphase situation. We note that the difference with the single- 
phase case occur only in Theorem 2.1, since finding the functions Bji„ui\{6,X) here is a more 
complicated problem. Let us return to the system (2.62), i.e. 

figj(X)S,[q](i)(0,X) + ^J,jjqj(0,X) = 0, (3.1) 

where the operator SE(X) and the functions A^,r i(0,X) are given by (2.63) and (2.64) respec- 
tively. It is easy to see that the operator B^^L (X) is a differential operator on the torus, which acts 
along a constant direction k°'d/d6°'. The corresponding vector field generates one-dimensional 
trajectories which are everywhere dense in the torus T™ for generic vectors k(X). For special 
values of k(X) the closures of the trajectories may be lower-dimensional tori, in particular, the 
closed one-dimensional trajectories in T'". Note that the operator B^,L{X) = B^,^,(U{X)) is 
exactly the Hamiltonian operator for the bracket (2.16) on the manifold of ?Ti--phase solutions A. 

The operator B^.iJJJ) has in general finite number of "regular" eigenvectors with zero eigen- 
values defined for all values of the parameters U and smoothly depending on the parameters. 
However, for special values of U the set of eigenvectors with zero eigenvalues is infinite and 
determined by the dimension of closures of the straight-line trajectories in T™, defined by the 
vector k(U).^ 

Let us define in the space of the parameters U the set Ai, such that for all U G A^ the tra- 
jectories of the vector field {k^(\J), . . . , /c™(U)) are everywhere dense in T™. From the condition 

rank ||9A;"/5C/''|| =m 

it follows that the set M is everywhere dense in the parameter space U and, moreover, has the 
full measure. 



^For some special brackets (2.16) the differential part can be absent in the operator _BE(U). The opera- 
tor _B^g, (U) reduces then to an ultralocal operator acting independently at every point of T™ . As a rule, the 
matrix B'!^, (U) is non-degenerate in this case. For example, for the ultralocal Poisson bracket 

{'ip{x),ip{y)} = iS{x-y) 
for the NLS equation 

iipt = ipxx + i^M^ip 

we have exactly this situation. Easy to see that system (3.1) represents a simple algebraic system in this case and 
is trivially solvable. The multiphase situation is not different here from the single-phase one, so all the arguments 
of the previous chapter can be used also for the multiphase case. However, for arbitrary brackets (2.16) the 
operators -B?p, (U) have more general form described above. 



36 A.Ya. Maltsev 

In the study of the solubihty of (3.1) we must first require the orthogonahty of the functions 
A^ir i(0,X) to the "regular" eigenvectors of SE (X) with zero eigenvalues. Let us prove here 
the following lemma. 

Lemma 3.1. Let K he a regular Hamiltonian family of m-phase solutions of (2.3) and (/ , . . . , 
I^) be a complete Hamiltonian set of commuting first integrals of (2.3) having the form (2.18). 
Let forlJ eM 

be the complete set of linearly independent eigenvectors of the operator i?E (U) on the torus with 
zero eigenvalues, smoothly depending on 6. Then 

1) the number of the vectors v^^, (6) is equal to the number of annihilators of the bracket (2.16) 
on the submanifold of m-phase solutions of (2.3); 

2) the functions yl^^,r , (0,X) are orthogonal to all the vectors ^\-[j(x)]i^)' ^•^• 

r>27r r2n 



/■ATT pZTT r\''^f) 



Proof. Consider the values of v^'JO) on any of the trajectories of the vector field k°'d/d6"' on 
the torus T™. According to the definition of regular Hamiltonian family of m-phase solutions 
of (2.3) the corresponding functions v^rjjJhx + Oq) should be the variational derivatives of some 
linear combination of annihilators of the bracket (2.16). We have then for a fixed value of 6q 



V 



«j(kx + 0o) = E4(U,^o)^ 



¥!=$(kx+0o,U) 

From relation (2.26) we have then 

v% (k(U)x + e^) = Y, «p(U, 0o) <(U) K Jij (k(U)x + 0o) . (3.2) 

By definition, for U € TW the trajectories of the field k°'d/d9" are everywhere dense in T™. 
Since both the left- and the right-hand parts of (3.2) are smooth functions on T™, we get then 
that they coincide on T™. We can put then ap(U, 6q) = ap(U) and write 

4u] (^) ^ E «p(U) n^(U) K^l^ (6) . (3.3) 

It is easy to see also that any linear combination of the form (3.3) gives a regular eigenvector 
of the operator B^,L(JJ(X)) with zero eigenvalue. 

Statement (2) follows then from relation (2.71) in view of representation (3.3). ■ 

However, despite the presence of Lemma 3.1, study of system (3.1) is much more complicated 
in the multiphase case if compared with the single- phase case. Thus, the presence of "resonances" 
for some values of the parameters U(X) may lead to insolubility of (3.1) in the space of periodic 
(in all 6°') functions. The set of the "resonance" values of U, as a rule, has measure zero. 
Nevertheless, the possibility of resonant parameters do not allow to transfer directly the methods 
of the previous chapter to the multiphase case. The following theorem shows, however, that the 
procedure of the bracket averaging is in fact insensitive to the appearance of the resonant values 
of U and can be used in most multiphase cases, as well as in the single-phase case. 
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Theorem 3.1. Let the system (2.3) be a local Hamiltonian system generated by the func- 
tional (2.17) in the local field-theoretic Hamiltonian structure (2.16). Let A be a regular Hamil- 
tonian family of m-phase solutions of (2.3) and {L^, . . . ,/^) be a complete Hamiltonian set of 
commuting integrals (2.18) for this family. 

Let the parameter space U of the family A have a dense setS C Ai on which the system (3.1), 
or, equivalently, the multiphase system (2.72), (2.73) is solvable in the space of smooth 27r- 
periodic in each 6" functions. Then the Dubrovin-Novikov bracket, obtained with the aid of the 
functionals (/^, . . . , / ), satisfies the Jacobi identity. 

Proof. As before, for a smooth compactly supported vector- valued functfon q(X) = {qi{X), 
. . . , qi\i{X)) we define the functional 

r'+OO 



/ + 00 
r{X)q,{X)dX. 
-oo 



Then, for arbitrary smooth, compactly supported in X and 27r-periodic in each 9" functions 
Q{0,X) = {Qi{0,X), . . .,Qn{6,X)) we define the functionals 

Qi{e,X) = Qi{0,X) 

-^,{o,v{x))M^''{ij{x)) ... Q^(e',x)^i,,{e',\j{x))——, 

Jo Jo y^'^) 

where the matrix M^'^(U) is the inverse of the matrix 

which is always defined according to the definition of a complete regular family of m-phase 
solutions of system (2.3). 

By definition, the functions Qi{0,X) are local functionals of U(A') 

Qi{e,x) = Qi{e,x,\j{x)) 

depending also on the arbitrary fixed functions Q(0, X). Everywhere below we will assume that 
Q(0,X) is a functional of this type defined with some function Q(0,X). 

Easy to see that the values of Qi{0,X) with arbitrary Q{6,X) represent for fixed values of 
the functionals U(Z) all possible smooth, compactly supported in X and 27r-periodic in each 0" 
functions with the only restriction 



r2n (\^^f) 

j^ Q,{e,X)%^{e,\^{X))j^-^=Q, VX,a = l,...,m. (3.4) 



t-2-K (.2tt 

Jo 
For the functionals Qi{6,X) we define the functionals 



^[Q] = C^ T- r a\o, X)Q. (^ + e, x) ^dX 



-oo JO JO 



{27ry 



Now, for fixed functions q(X), p{X), and functional Cl{d,X) consider the Jacobi identity of 
the form 

{g[Q] , { -/[q] , >^[p] } } + { -^^Ip] , im] ' ■^[q] } } + { -^[q] ' { -^[p] ' m] } } = 0- (3-5) 
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Expanding the values of the brackets {J'^{X), J^'(y)} in the neighborhood of the submani- 
fold /C, as in the single-phase case, we can write 



{r{x),nY)} = {r{x),j^{Y)}\^ 

r /'27r /'27r 

Jo 



s>0 



6ip'^{e,w)Jo 



A7{^{e',X),e^x{0',X), 



dre' 

(2^ 



K 



dre 



X 9\e, W)j^-^dWe^5^^\x -Y) + 0(g2), 

where all the values on the submanifold K, are calculated at the same values of the functionals 
[U(Z)] as for the original function (p(6,X). 
Let us introduce by definition 



^{J[ci], J[p]} 



K 






s>Q ■ 



(27r)" 



(3.6) 
2u{X)p^^sX'^X. 



Note that the notations 5/5g {0,W), generally speaking, are not natural in our situation 
because of the dependence of the chosen system of constraints. Nevertheless, the preservation 
of these notations can better clarify the algebraic structure of the further calculations. 

The values defined by (3.6) can be represented in the form of the graded decompositions at 
e — 7- on the submanifold /C. By virtue of (2.35) it is easy to conclude that the expansion in e 
of the quantities (3.6) begins with the first degree in e 



<5{J[q],J[p]} 



Sg'^ie^W) 



'5{J[q],J[p]} 



K 



5g\e,W) 



+ e' 



'^{«^[q],'^[p]} 



/cm 



5g\e,W) 



+ ••• . 



/C[2] 



(3.7) 



The leading term of (3.7) can be divided into two parts, corresponding to the functions 
eQ^((^, ey^x, . . . ) and A\^{ip^e(fx-,---) and containing the quantities quiW)p^{W) and 
qv{W)Pf_i,w{W) as local factors, respectively. 

The values of { J[qi , Jrpi } , are obviously invariant under transformations of the form 



V?(0,X) ^(^(6' + A0,X). 



(3.8) 



From the invariance of the functionals U(X) in such transformations, we can write for the 
corresponding increments of the constraints (2.32) 

5g^{e, X) = /(6> + 56, X) - /(6>, X). 

As a consequence, we can write on the submanifold /C 



27r 



,0 JO Sg\e,W) 

a = 1, . . . , m. 



$ 



e J (27r)"^ 



(3.9) 



The above relation is satisfied to all orders in e. By the arbitrariness of the functions qy{X), 
the relation (3.9) in the leading order can be strengthened. Namely, according to the remark 
about the form of the leading term of (3.7) we can write for any W 



Jo Jo 



Sg''{e,W) 



IC[1] 



^U^-^ + sMw? " 



(27r)" 



0, 



a 



,m. (3.10) 
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Quite similarly, we have the relation 



+ / 5'(0,X)Q,,,. [^ + e,X^ lsgn(X-y){A;"(y),J[q]}^^dydX 






S(X) 



l>0 

. fS{X) 



d'^e 



, +^Mx)){u^x),j^^]}j^-^dx 



JQ.{ 



+ e,x\^lJ^^ + e,v{x)\ 



1 (V^fl 

X - sgn(X - Y) {k-{Y\ J[^]} 7,^— ^dydX 



2e 



(2vr)' 



+ / g\e,X)Q,^e^ ( ^ + 0,x) lsgn(X - y) {A:"(y), J[q]} ^^dFdX 



+ 



S(X) „ „\ ,„„.„. . , d™0 



jg\e,X)Q,^u^ (^^ + 0,Xj {t/-(X),J[,]}^dX 



According to relations (3.4), we can actually see here that the third term in the above 
expression is identically equal to zero. 

According to the form of the constraints, we have again in the neighborhood of /C 



{m],J[q\\ = {5[Q],-/[q]}|yc 



+ 



Qi 



S(X) 



+ e\xW, 



6C}^J^{e',X),e^x{0',X),...) ^mff 



r^{i) 



l>0 



6^''{e,w) 



(27r) 



q^,lx{X)dX 



Q,{e',xWu.{e\\j{x)) 



d^e' <5{t/-(X),J[q]} 



(27r)"^ 8ip^{e,W) 



dXxg\e,W)--—dW 
K (2^) 



+ / QkM<^ 



+ / Qi,^ 



S{W) 



)1 (Y'^f) 

-^g^iW -Y) {k^{Y),J^^^]\^dY X g\e,W)j-^dW 

)ATnQ 



provided that all the values on the submanifold /C are calculated at the same values of the 
functionals [U(Z)]. 

For the quantities {^[q], <^[q]} we can introduce by definition 



5g^{e,W) 



= ^Q^\^^ + e\x 



E 

«>0 



^5C'^^{<^{e\X),e^x{9\X),. 






-%,ix{X) 



K 



d^e' 

(27r)"^ 



dX 



dX 



K. 
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The quantities S{g[Q], J\f^}/5g^{0, W)\k: have the order 0(1) at e ^- 0. We have also 

c'>^{^{e',x),e^x{e'.x),...)^s'^^{^{e\x),e^x{e',x),...) 

according to relation (2.37). 
Let us introduce the functions 



(3.11) 



S'^^^{ip,ip^,...\ 



95^'^(v,<^x,. 



dip'' 



Ix 



Using now Lemma 2.3 and relations (3.7) we can write for the leading term of (3.11) 



^{9[Q\,Jh\ 



5g^{e,W) 



m 



Y{-l)^k^^{W) ■ ■ ■ k'^'iW)——^^—— 

ft ' ^'^''' 



l>0 



+ e,w] sf^^^^{^{e,w),e^w{e,w),...\ 



KV 



+ io''>'{W)q^{W)Qk,e'^ ( ^^ + e,W ,, 



(3.12) 



where 



sf,.A<^{e,w),e^w{e,w), 



'W) 



K\i 



Note one more property of the values JjgijQ], J[q]}/55^(0, VF)|y(;[o]. As we saw earlier, the 
values {(7[Q], -/[q]} are of the order 0(e) at e — ;■ on the submanifold /C. This property is 
preserved also under the overall shift of the initial phase (3.8). 

Indeed, for 



"P 



e, X) = ¥ (^^ + e + Ae, v{x)] 



(3.13) 



we can write 



{5[Q],«/[q]} - J Q^ 



^^ + e,z)^sgn{Z-W) 
e / 2e 



{k'{W). J[,| }<!.'( ?^ + fl + AS, 0(2)) ^dZdW 



+/fti^+«,z)car*^^*"' 



+ 9 + Ae, V{Z)\ , ^$0/^+ 6» + A6>, V{Z)\ ,...") 



>'Qf^iZ)—^dZ + 0{e). 



(2vr) 

Because of the invariance under translations (3.8) the value of the bracket {k°'{W), -/[q]} on 
the functions (3.13) is equal to its value on /C 



{k-{W),J[^]} =e{u-''iW)q4W))^ + 0{e'). 



Whitham's Method and Dubrovin-Novikov Bracket 



41 



Similarly, we have on the functions (3.13) 



cz ( * ( ^ + 9 + Ae, v{z)\ , F*e. (^ + e + Ae, u(z) ) . . 



^(0) 



We then obtain 



S(Z) 



+ 6> + A6',U(Z) . 



/O 1 777, /I 

<?^(Z)a;"^(Z)— [Q,ie, Z)^\e + A0, U(Z))] ^2^dZ + 0(e) = 0(e). 



As a consequence, we can write 
, ra 



'0 

a = 1, 



$i ( ^^ + 0, U(VF)1 T^^^diy = 0, 



(2vr)^ 



(3.14) 



for the main part of (5{(7[q], J\^lbg^{Q^ VF) on /C. 

Using again the fact that the relation (3.14) contains arbitrary functions q^i{W\ appearing 
in the integrand expression in the form of local factors, we can rewrite (3.14) in a stronger form 



2-K 



^" ^{g[Q],>/[q]} 

^g\e.w) 



^ 



^M + ,,u(M.)^ d™^ 



(2^) 



0, 



(3.15) 



VM^, a = 1, . . . , m. 



We now turn back to the Jacobi identity (3.5) for the functionals yrqi, J[q], and Jrpi. It is not 
difficult to see that after the restriction on K the leading term (~e) of the relation (3.5) can be 
written as 



/{'wi.^'C'"')}!..,!^ 



d"^6i 



K{\\ 



(2vr)^ 



AW 






{J[q],/(0,t^)}|^ 



'^{5[Q],^[p]} 



W (^5'=(6',Ty) 



d"g 

yc[0](2vr)- 
d'"6' 



K\i 



{27ry 



dW 



-dW = 0. 



The above identity can again be written in a stronger form. Namely, making the change 
Qi{9,W) -^ Qi{e,W)i2i{W), we get the corresponding change, Qi{e,W) -> Qi{e,W)fii{W), 
where ^J-iiW) are arbitrary smooth functions of W. By virtue of (3.4), (3.10), and (3.12) it is 
easy to see then that the integrands are smooth functions of and W, containing fJ-i{W) in the 
form of local factors. By the arbitrariness of fii{W), we can omit the integration over W in the 
above integrals and write for every W 



2-K 



+ 



27r 



27r 



/'27r f-l 

Jo Jo 



|5[Q],5 (0, W)||^[o] 5g>^ie,W) ^[1] Wr 



,l"^^M'''^'''^^^'^W ^,.(^,Ty) 



^[0](2^ 
d™^ 



/C[0] 



(2vr)' 



0. 
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Finally, using the relations (2.61) and (3.15), we can write the above identity as 



f2-7r /•2-7r 



/ •••/ {5[Q],/(0,W^)}| 
Jo JO 



'^{^[q],>/[p]} 



'CM 5g''{e,W) 



2lT f2-K 



Jo 

1-K /-StT 



4*(^.<'.«' 



+ / ■■■/ -^^ 

JO 



[i][q] 



e 

S(T^) 



d™0 

yc[i](2^ 



+ 6',iy 



55^^(0, W^) 



bg^{e,W) 



d"6' 
d"6' 



ycfoi (2vr)^ 



0, 



where the functions ^f^^,r , (0, W^ are introduced by formula (2.64). 

Now assume that the values A^^^,^(^^W^ for U(H^) G 5, according to (3.1), can be repre- 



sented in the form 



ni][q]^ 



with some smooth in 0, 27r-periodic in each 0° functions B j\rA(\\{0 ^W\ 
We can then write for U(Ty) G 5 



(3.16) 



27r r2-K 

JO 

+ 






K.{\\ 



{27ry 



(3.17) 



Jo 

2tt f2iT 



H9[Q],Jlci]} 



^0 



^roifsiW^.Mi)(^^+^'^ 






ycfoi(2^ 



where 



S(VF) 



KH 



d™^ 

(27r)™ 



= 0, 



^io][s] (W) = E i3g) (* (^ + 0, U(T^)) , k^W)^e^ (^ + 0, U(H^)) , . . . ) 



s>0 



QS 



09'^! ■ ■ ■ 06""' 



Using the expression (2.54) for the bracket of constraints on /C, as well as the relations (2.42), 
(3.4), (3.10), we can obtain also the following relation 



2tt i-2n 

Jo 



{5[Q],<7 (0,W^)|l^[o] Sg''{0,W) ^[1] (2^ 



Jo 



[0][S]* 



^{Jlci],J[p]} 



6g'^{e,W) 



m (2^' 



Expanding also the remaining terms of the identity (3.17) according to (3.12), we get for 
IJ{W) G S 



2tt r2n 



Jo 



B^^,^,iW)Q,{^^ + e,W 



^^ms]' 



^{«^[q].^[p]} 



27r r2iT 



d'"6' 



^M[s](^)^.-[pl( 



1) 



6g>'{e,W) 

S{W) 



Jo (2vr) 

Yi-lYk^'HW) • • • k'^'iW)--^^-—— 






l>0 
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+ io''^'{W)q^{W)Qk,e'^ ( ^^+e,W 



K\^ 



Q,iW) 



+ 



2tt r'2-iT A^^O 

"Jo (2^ 



^asiW^.M(i)(^^+^'^ 



Ql 



Yi-lYk'^^iW) ■ ■ ■ k^'iW)--^^—— 



+ e,w]sl^^,^{^{e,w),...) 



+ u;''''{W)p^iW)Qk,e'^ { ^^ + 0,W 



K\i 



l>0 

0. 
Consider now the Jacobi identity of the form 

{9lF],{g[Q],Jlci]}} + {9lQ],{Jlci],9lF]}} + {Jlci],{g[F],9lQ]}}=0 



(3.18) 



(3.19) 



for arbitrary fixed functions q(-^), and the functionals P{6,X) and Q{6,X) defined as before 
with the aid of arbitrary functions P{6,X), Q{6,X). 
According to the relations 

{9[Q,P],U^{W)}\^ = 0{e), {j[q],/(0,VF)}|^ = O(e), { J[q], [/^(TV)} |^ = 0(e), 

it's not difficult to see that after the restriction on /C the major term (in e) of (3.19) will be 
written as 



{5[P],/(0,Vr)}| 



H9[Q],J[ci]} 



'CM 6g''{e,W) 



d'^e 



(2vr)' 



-dW 



{9lQ],9H0,W)} 



H9[P],Jlq]} 



ICli 



dgi^ie^w) 



d'^e 

(27r)™ 



dT^^O. 



Again recalling that q^iW) are arbitrary functions of W appearing in the integrand in the 
form of local factors, we can write the above relation in a stronger form. That is, for every W 



" '"' y„,,\o,w)}\ ''"M'-'w' 



JO 



d'^e 



27r r-2iT 

^0 



'CM 6g''{e,W) 
{9[Q],9 [d,W)}\ 



6g''{e,W) 



jcm (2^)' 



0. 



(3.20) 



As well as in the case of the identity (3.18), using the relations (2.54), (2.42), (3.4), (3.15), 
we can write the identity (3.20) in the form 



2tt r2'K AiTT-f) 

"Jo WT 

X 



m 



B^miw)Pj 



S{W) 



+ e,w 



Ql 



Yi-lYk'^^iW) ■ ■ ■ k'^'iW)--^^—— 



l>0 



+ e,w]sl^^,^{^{e,w),...) 



IC\i 



q>^{w) 
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+ io''^'{W)qf,{W)Qk,eo. ( ^^+e,W 



2tt f'2iT f^ni/D 

"Jo Wr 



jkj 



S{W) 



Blfl^[s]iw)Qj[^r^ + e,w 



l>0 



^ + e,w^sl^^,^i^ie,w),...: 
+ uj'"'{w)qf,{w)Pk,e'^ (^^ + e,w 



/C[0], 



0. 



(3.21) 



Note now that the values of Q(0,X) and P{0,X) are arbitrary 27r-periodic functions of 9, 
satisfying the conditions (3.4). In particular, we can put in (3.21) at U(VF) G S 

P(0,T^) = B[p](i)(0,VF) or P(0,W^) = B[q](i)(0,t^). (3.22) 

By analogy with (3.12) we introduce for convenience the notation for U(W^) G 5 

'5{5[eB[p,(i)],'/[q]} 



5g\e,W) 



/C[l] 



Y{-iyk^^{W) ■ ■ ■ k'^'iW)—-^—— 



%{W) 



+ a;-'^(W^)g^(t^)i?fc[p](i),,. [^^ + e,W 



(3.23) 



for arbitrary smooth functions q(-'^), p(^)- 

Note that the functional g[eB\ in,] is not defined on the whole functional space, so the rela- 
tion (3.23) plays just a role of a formal notation for U(VF) G 5. 

Using now (3.21) for the functions (3.22), we can rewrite (3.18) in the form 



/•27r j- 

Jo Jo 



2-K jm 



(TO 



jkj 



/o Jo (2vr) 

HJ[ci],J[p]} 



BomiW)Q^ 



S{W) 



6g''{e,W) 



provided that U(VF) G S. 



[0][S]. . -. , g 



+ e,w 



JC[1] 



SgHO,W) 



+ 



<5{ff[.B[q,(i)],-f[p]} 



K[l] 



6gH0,W) 



IC[1], 



Using the skew-symmetry of the Hamiltonian operator BJ-,r„JW) on A we can then write 



[0][S]^ 



/•27r p 

Jo Jo 



Jo (2vr) 



-d-.^7S(H0^,^^ 



^S[S](^) 



'^{^[q],'^[p]} 



dgf'ie^w) 



'^{5[eB,pjn)],'/[q]} 



K[l] 



5g''{e,W) 



+ 



'^{5[.B[„](i)],-^[p]} 



JC[1] 



SgHO,W) 



K[l]. 



0. 



The values Qj{0,W) are arbitrary smooth 27r-periodic functions of 6 with the only restric- 
tion (3.4). We know also that the values in the brackets are smooth 27r-periodic in each 6*" 
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functions of 6 for U(H^) £ S. As a consequence, we can write for U(VF) G S 






5{%Brpin)],-f[q]} 



1C[1] 



dgi^ie^W) 



+ 



Hff[6B[qj(i)],^[p]} 



^[1] 



5g^{e,W) 



K.[l\, 



= Y,o^''{^{W),\]w{W)We^ 



a=l 



S{W) 



e,v{w) 



with some coefficients a"(U(H^),Uiy(VF)). 

The values in parentheses are smooth 27r-periodic functions of at U(VF) G 5. At the 
same time the trajectories of the vector field (A:^(VK), . . . , /c'"(W^)) are completely irrational 
windings of the torus T*". We can therefore say in the case of a regular Hamiltonian family A, 
that up to a linear combination of the regular eigenvectors v^^.^s, (S(VF)/e + 0) oi BL,,^AW) = 

-Bj^Qirg, (U(VF)), corresponding to zero eigenvalues, the value in parentheses is a linear combination 
of the variational derivatives (2.7), generating linear shifts of the phases on A. For complete 
Hamiltonian set of the integrals (I^, . . . ,1^) we can then write according to (2.8) and (3.3) 



^{J[ci]^ J[p]} 



5g^{e,W) 



^{5[eB[p]n)],^[q]} 



1C[1] 



5g>^{e,W) 



+ 



'^{5[eB,qjn)],^[p]} 



1C[1] 



5g^[e,W) 



K[l] 



= 5^5,(U(T^),Uh.(T^))4V)] i^+o) 



(3.24) 



with some coefficients hq{\]{W),\]w{W)) at \]{W) £ S. 
Consider the Jacobi identity of the form 

{{J[ci], J[p]} , J[r]} + {{J[p]^j[r]} , ^[q]} + { { -^[r] , -/[q] } , ^[p] } = 

with arbitrary smooth functions q(A), p(A), r(A). 

In the main (~e'^) order on /C given identity leads to the relations 



^U[cih Jlp]} 



5W{W) 



JC[1] 



{U^W),J[r]}\i^,,.dW + C.p 



+ 



II 



2-K I-2TX 



^{Jh^J{p\} 



bg^(e,W) 



/C[l] 



{/(0,t^),J[r]} 



d™^ 

^\^\ (27r)™ 



dl^ + c.p. = 0. (3.25) 



Again, using the relations (2.61), (2.64) and (3.10), we can replace the identity (3.25) by the 
following relation 



^{Jh^J\p^ 



5f/^(VF) 



/C[l] 



{t/^(VF),J[r]}|^,ydTy + C.p 



+ 



2-K I-2TX 







^{Jh^J{p\} 







bg^(e,W) 



/C[l] 



^[l][r] 



' 27rV» ^ 



Using relations (2.71) and representations (3.16) and (3.24), we can write for U(M^) G S 
io 70 <5ffH^,W^) 



A 



/C[l] 



[l][r] 



S(PF) 



+ 6>,VF 



d™0 

(27r)™ 



+ c.p. 



2-K (■2-K 

^0 



^[o][s]W^.Mi)^^^ + ^'^ 
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'^{5[.B[q](i)],>/[p]} 



5g''{e,W) 



H^BrpKnl'-^Wi 



JC[1] 



6g''{0,w) 



K[ll 



+ 



2-K 



2-K 



^fo][S](^)^i[q]{l) 



S(VF) 



e,w 



'5{5[.B[pj(i)],^[r]} 



bg^(Q,W) 



H5[.Bm(i)],>/[p]} 



/C[l] 



5g'=(0,VF) 



K.\\\ 



+ 



27r 



27r 



,^. mno__ /'S(t^) 



^folfsi(^)^i[p](i) 



+ 6>,VF 



'^{5[.B[,](i)],^[q]} 



bg^(e,W) 



'5{5[eB[q,(i)],'/[r]} 



^[1] 



5g'=(0,H^) 



/C[l]. 



(27r)" 



(27r)™ 



(27r)™ 



(3.26) 



Similarly to earlier arguments, substituting now in the identity (3.21) 

q.{e,w) = B[r](i)(0,T^), Pie,w) = B[p](i)(0,w^) 

for U(W^) G 5 we obtain the identity 



2n 

^0 

2tt 



2-K 



^[01[S]W^.[P1(1)(^^+^'^ 



•^I^BirKDl'-^Iq]} 



/Cfl" 







27r 



^fo][s]W^.Mi)f^® + ^'^ 



'^{5[6B[pj(i)],^[q]} 



(27r)™ 



55'=(0,iy) 



/C[l] 



(27r)™ 



(3.27) 



0. 



Using the cyclic permutations of the functions q(-'^), p(^), and r{X) in the identity (3.27), 
it's not difficult to see that the right-hand part of the relation (3.26) is identically equal to zero 
at U(PF) E 5. It's not difficult to see also that the left-hand side of the expression (3.26) is 
a smooth regular function of the parameters U(VF). Using the fact that the set 5 is everywhere 
dense in the parameter space U, we can conclude that the left-hand side of equation (3.26) is 
identically equal to zero under the conditions of the theorem. 

We have, therefore, that under the conditions of the theorem, the identity (3.25) implies the 
relation 



5W{W) 



K.[l\ 



{[/^(W^),J[r]}|^[ydVF + C.p. = 0. 



Using the relations 



^{J[ci], J[p]} 



5W{W) 



K[l\ 



5U^{W) ' 



{^^(W^),J[r]}|^j,, = {C/^(Ty),C/[r]} 



DN' 



we obtain the Jacobi identity for the Dubrovin-Novikov bracket on the space of the func- 
tions U(X). ■ 

Remark 3.1. According to relations (2.47) the functionals 

f+OO 



i^° 



A;"(U(X))dX 



are annihilators of the Dubrovin-Novikov bracket (2.21) since we have {if",C/rqi} = for any 
functional t/[q]. As a consequence, we can claim that the functions k°'{\J) represent a part of the 
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flat coordinates for the metric g'^^(\J) = {A'^^) connected with the Dubrovin-Novikov bracket 
for the Whitham system. Besides that, we have 

{k'^{X),k^iY)}^^ = 

according to relations (2.48). 

The remaining part of the flat coordinates of ^'^'^(U) and the corresponding annihilators of 
the bracket (2.21) are defined by concrete form of the initial bracket (2.16) and the family of 
m-phase solutions of system (2.3). 

Let us prove now the second theorem justifying the invariance of the Dubrovin-Novikov 
procedure. 

Theorem 3.2. Let A be a regular Hamiltonian family of m-phase solutions of (2.3). Let 
{I^, . . . ,L^) and (/'^, . . . , /'^) be two different complete Hamiltonian sets of commuting first 
integrals of (2.3) having the form (2.18). Then the Dubrovin-Novikov brackets obtained using 
the sets {I^, . . . , L^) and [L'^, . . . , L'^) coincide with each other. 

Proof. Note that the sets (/^, . . . , I^), {I'^, ■ ■ ■ , I'^) correspond to the two systems of coordi- 
nates (C/^, . . . , U ), {U'^ , . . . , U ) on the family A, given by the averages of the functionals I 
and I'. We have to show then that the Dubrovin-Novikov brackets obtained using the set 
(I'"*^, . . . ,/ ) coincides with the bracket, obtained using the set (I"*^, . . . ,/ ), after the corre- 
sponding change of coordinates, 

U'" = C/"^(U). 

Consider the values of the functionals 

and the values of constraints g^{0,X), introduced with the aid of the functionals J{X) by 
formula (2.32), as a coordinate system in the neighborhood of the submanifold /C. 

For values of the functionals J'^{X) on the submanifold /C we then have the relations 

J-(X)|^ = [/- (J(X)) + Y.e'j'(;^ (J, Jx, . . . ) , 

where j'X-. are smooth functions of (J, Jx, • • • )) polynomial in the derivatives and having degree I. 
Expanding the values of J"^{X) in the neighborhood of the submanifold /C, we can write 



J-(X) = [/-(J(X)) + ^e'i'(r)(J,Jx,...) 

/+00 /■27r /•2-iT rl"^/9 

-oo I "Jo ^"(^'^'^'^)5^(^'^)(^^dy + o(g2), 

where, according to the form of the constraints (2.32), we can put 

Tr{X,e,Y,e) = Y,K(i)MO,X),e^xie,X),...) e^6(^\X-Y) 



l>0 



^^{^{e,x),e^x{o,x),...: 

i>o '^i^ 



e'5«(X-y). 
K 
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Considering the functionals J', , = J qu{X)J"^{X)dX with arbitrary smooth (compactly sup- 
ported) functions qi^{X), we can write in the vicinity of /C 






i>i 



l>0 



'A. 
dy' 



dX 



q.{Y)Il'^,^{^{e,Y),e<fyi9,Y),.. 



K. 



(3.28) 






The leading term (in e) in the second part of the expression (3.28) is given by the expression 



jq^iY) Y^[-l)^k'^^{Y) . . . k"'iY)U\1,^,^,...,., U 



l>0 

and coincides with the value 

Jo 

where 

r f2iT f2n 



m,e.r)..^y,e.,^^.r 



a: 






AITIQ 

+ e]g\e,Y)7^-^^Y, 



<S(^) 



V(^) 



Jo Jo 



d™6i 

(27r)™ 



(^(6l)=*(6»,U) 



As the values Qfui(^)' ^^^ values QfuU^) represent regular left eigenvectors of the opera- 
tor i^fui corresponding to the zero eigenvalues. In the case of a complete regular family of 
m-phase solutions, we have therefore 



c:[u5w=Er.^(u)^%w 



for some functions rj^(U). 

We can write, therefore, up to quadratic terms in g{9, X) 

j q,{X)J'\X)dX = jqAX) U'^ (J(X)) + ^ e'i|r) (J, Jx, 



dX 



+ 



Jqu{X) 



Er^(uK?V^i(^^ + ^)+o(.) 






Consider the Poisson brackets 



{-^[ql'^llL = [j qAX)r{X)dX, j p,{Y)J'^{Y)dY 
= I qu{X)^{X)\^S^j\X), I p,{Y)J'^{Y)dY 

+ y'g,(x)j;r^^(u(x)) 



/c 



K. 



+ 0{e') 



dX 



M 



S(X) 



[U(X)] I ^ + ^ j + O(^) 

{5X0,x),j'^(y)}|^p,(y)^^-^dXdy. 



Whitham's Method and Dubrovin-Novikov Bracket 49 

By Lemma 2.4' we have the relation {g^{0,X), J! A\x: = 0(e) on the submanifold /C. In 
addition, completely analogous to the relation (2.71) holds the relation 



f 

Jo 



l"''^mi{^ + '>)U-(o.^UU}lm^.^» 



r2TV /■27r 

/O 

by virtue of the original dependence of the constraints g^{6,X). We thus obtain 



.(X)|^(X) \ J\X), I p,iY)J'^iY)dY 



K. 



jj q,{X){r{X),j'^^{Y)}\^p,{Y)dXAY 

dX + 0{e^). 

Repeating the arguments for the functional J p^{Y)J'^^{Y)dY we finally obtain 

jj q,{X){j"'{X),J'^^{Y)]\^p,{Y)dXdY 

,{X)^^{X){J\X), J7(y)} |^__(y)p^(y)dXdy + 0{e'). (3.29) 

Given that the principal (in e) terms in the expressions {J'^{X), J''(Y)}\fc and 
{J"^{X), J'^{Y)}\k: coincide with the Dubrovin-Novikov brackets, obtained with the aid of the 
sets {I^, . . . , /^) and {I'^, . . . , I'^) respectively, we conclude from (3.29) that 

{[/-(X),[/'^(y)};^ = ^^{X){U\X),U^{Y)}^^^^iY), 

which means the coinciding of the brackets {•, -Jdn and {•, -I'dn- ■ 

Finally, we prove the theorem about the Hamiltonian properties of the Whitham system (2.10) 
under the same conditions as before. 

Theorem 3.3. Let A be a regular Hamiltonian family of m-phase solutions of (2.3). Let 
{I^, . . . ,L^) be a complete Hamiltonian set of commuting first integrals of (2.3) having the 
form (2.18) and H be the Hamiltonian function of the system (2.3) having the form (2.17). 
Then the Whitham system (2.10) is Hamiltonian with respect to the corresponding Dubrovin- 
Novikov bracket (2.21) with the Hamiltonian function (2.22) 

f+OO 



/+0O 
{Ph) (U(X)) dX 
-oo 



Proof. By Theorem 3.2, without loss of generality we can assume that the Hamiltonian func- 
tional H belongs to the set (/^, . . . ,/ ), H = I^° . It is easy to verify that the corresponding 
Hamiltonian H"""" generates in this case the system 

u;^ = ^[{Q''^'){^{x))], 

i.e. exactly system (2.10). ■ 

As the simplest example of Theorems 3.1-3.3 consider the procedure of averaging of the 
Gardner-Zakharov-Faddeev bracket for the KdV equation. 
Consider the Gardner-Zakharov-Faddeev bracket 

Mx),ip{y)} = 5'{x-y) (3.30) 
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for the KdV equation 

^t = V^x - ^xxx 

with the Hamiltonian functional 

As is weU known, the KdV equation has a family of m-phase solutions for any m, given 
by the Novikov potentials, which represent the stationary points for the higher KdV flows. 
According to the Hamiltonian structure (3.30) it is also equivalent to the extremality of a linear 
combination of higher integrals of the KdV 

ciJ/i + C25/2 + • • • + C„+2W™+^ = 0, (3.31) 

where 

I^ = N = f ifdx 



I^ = P= / ^dx 



is the annihilator of the bracket (3.30), 

2 

is the momentum functional of the bracket (3.30), I^ = H, and I^ , k> A are higher integrals of 
the KdV equation. 

As was shown in [41], all the systems (3.31) are completely integrable finite-dimensional 
systems with quasiperiodic solutions in the general case. As is well known, the correspond- 
ing m-phase solutions of KdV can be represented in the form (1.2) with the theta-functional 
expressions for the functions $(0,U), satisfying the system [10-13,17,26,27] 

w°(U)$0« = A;"(U)$$ec - A;"(U)A;^(U)F(U)$0„e^0^. (3.32) 

The parameters of the solution (i?i, ..., £'2^+1)1 Ei < E2 < ■ ■ ■ < -E2m+i represent the 
branching points of the Riemann surface of genus m, and, together with the initial phases 
(0Q, . . . , 0^), completely determine the corresponding solution ip{x^ t) [41]. 

The eigenmodes of the linearized operator (3.32), as well as the adjoint operator, were studied 
in detail [4,5,28-30]. In particular, we can state that the families of ?TT,-phase solutions of KdV 
are complete regular families in the sense of Definition 1.1 and the corresponding systems (2.10) 
for any independent set of 2m + 1 higher integrals of KdV (/^ , • • ■ , 1^™-+! ) represent the regular 
Whitham systems for the KdV equation in the sense defined above. It is well known [22,49], 
that the Whitham system for KdV in ?7z-phase case can be written in the diagonal form where 
the parameters -Ei, . . . , E2m+i are its Riemann invariants. 

Regular zero modes of the linearized operator (3.32) on the torus T™ are given by the functions 
$0a(0,U), 0" = 1, . . . ,m, and by the function d^/dn (at constant k and u), where n is the 
value of the functional N on the corresponding solution. Similarly, regular zero modes of the 
adjoint linear operator on the torus T™ are given by any m + 1 linearly independent variational 
derivatives of the higher integrals of KdV on the family of the functions $(0, U). 

The families of Novikov potentials are obviously regular Hamiltonian families with respect to 
the bracket (3.30), while any set of 2m + 1 independent higher integrals of KdV is a complete set 
of commuting functionals (2.18) on such a family. We also note that the variational derivative 
of the annihilator of the bracket (3.30) always appears on the family of ?TT,-phase solutions in the 
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form of a linear combination of ttt. + 1 variational derivatives of higher integrals of KdV according 
to the original construction of the Novikov potentials. 

Investigation of the system (3.1) can be carried out in this case as follows. We first note that 
the operator Biqi (X) is given in this case as 

and, thus, is the derivative along the vector field {k^{X), . . . , k'^{X)) on the torus T"*. 

The regular zero mode of the operator Biq-]{X) is a constant function on the torus and is 
orthogonal to the functions ^[i][q] {^i X) by Lemma 3.1. This fact is easily verified also by direct 
computation. Indeed, the functions Ami^UOjX) have in this case the form 



Al][q](^+^'^)={¥'(^.^)"^[q]}L[l]-^^''(^ + ^'U(X)){^'^(X),J[q]}|^ 



[!]■ 



Let us assume for simplicity that the functional 

N{X)= ... ^(0,X) ^ 

Jo Jo 



{27Ty 



is included in the coordinate system \J{X), N{X) = U^{X), which, by Theorem 3.2, does not 
affect the Dubrovin-Novikov bracket. We have then 






It is easy to see then that 



2-K i>2t: r\"^f) 



JO 



^[l][q](^, 



The operator B\q\{X), however, has also irregular zero modes which arise in the cases when 
the closure of trajectories of the vector field {k^{X), . . . , kP^^X)) on T*" represent the lower- 
dimensional tori T^ C T™. Modes of this type arise obviously on a set of measure zero in the 
parameter space U, i.e. \\J}/Ai. 

We assume here, without proof, following easily verifiable from the theta-functional represen- 
tation fact. The Fourier harmonics ^^^["^""(X) of the functions A\i^\^{6,X) in the variables 9 

Ai][q](^'^)= E A^^-{X)e^v{^nie^ + ■■■ + ' 

ni,...,nm 

decay faster than any power of |n| at |n| — )■ oo, where 



trimOn 



i2 I . . . I Ti2 



|n| = A/nf H hn^. 



Let us define the Diophantine conditions for an arbitrary set of values (^^(U), . . . , /c'"(U)). 
The vector (fc^(U), . . . , /c'"(U)) is a Diophantine vector with the index u > and the coefficient 
A > 0, if 

|niA;i(U) + --- + n^A;'"(U)| > yl|n|-'^ 

for ah (m, . . . , n^) G Z- ((m, . . . , n^) / (0, . . . , 0)). 

Let us use now the following well known theorem (see, e.g., [3,45]): 
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For u > m — 1 the measure of the set of non-Diophantine vectors {k^, . . . , /c™) in M™" is equal 
to zero. 

Using the condition rank ||(9/c"/5C/^|| = m we can also state that the measure of the corre- 
sponding parameters U, such that k(U) are non-Diophantine vectors with index v > m — 1, is 
also zero in the parameter space U. 

We can now define the sets S^ in the space of the parameters U, such that the wave vectors 
(fc^(U), . . . , A;'"(U)) are Diophantine with index z/ if U G S^. In this case all 3,^ with index 
ly > m — 1 are everywhere dense, and also have the full measure in the parameter space U. 

It is easy to see also that the system (3.1) is resolvable in the space of 27r-periodic in all 
9" functions for each JJ & S^- Indeed, due to the absence of the zero Fourier harmonic of the 
function A^ij jq] {6, X) and conditions of decreasing of Fourier coefficients ^mi„r'" {X) we can put 

Tjni...nm.( Y\ _ Ani...nm( Y\ 

"^WW ^^^"mifci(U) + --- + m^fe-(U) WW ^ ' 
for (ni, . . . , Tim) / (0, . . . , 0) and restore the solution 

ni,...,nm. 

(ni,...,nm) / (0, ...,0) 

as a smooth m-periodic function of 6. 

Assuming, therefore, S = S^ for any v > m — 1, we meet all the conditions of Theorem 3.1. 
We can thus formulate the following proposition. 

Proposition 3.1. The Dubrovin-Novikov procedure is well justified in the averaging of the 
Gardner-Zakharov-Faddeev bracket on the m-phase solutions of KdV for any m and provides 
a local Hamiltonian structure of hydrodynamic type for the corresponding regular Whitham sys- 
tem. 

We note here that both local and weakly nonlocal Hamiltonian structures of the Whitham 
hierarchy for KdV were investigated in detail in the papers [2,44]. 

In the more general case the constructing of the set S requires the study of the eigenfunctions 
and the eigenvalues of the Hamilton operator B]L{\]) on the manifold of Tn-phase solutions 
depending on the values of the parameters U. For m-phase solutions, given by the algebraic- 
geometric families, in this case are very convenient methods similar to those used in [4,5,28,30,47] 
to study the spectrum of auxiliary linear operators of integrable hierarchies. 
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